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It has long been known that pipe bends are much more flexible than 
can be explained by the ordinary formula for bending of curved bars. 
Experiments by Professor A. Bantlin of Stuttgart? show that the displace- 
ment of the ends of double-offset pipe-bends, when subject to forces tend- 
ing to alter the curvature, may be as much as five times as great as found 
by the ordinary theory. Professor Bantlin ascribed this to corrugations 
on the concave side and stretching on the convex side produced in man- 
ufacture, but Dr. Th. v. Karman of Géttingen advanced the theory that 
the principal cause of the discrepancy was a deformation of the pipe by 
which the tube section under bending takes a slightly flattened, oval 
form. 

It can be seen at once that when a pipe is subject to a bending couple, 
tending for instance to increase the curvature, the stress forces will have 
a tendency to push the parts on both sides of a transverse section towards 
the neutral surface. Hereby the material farthest from the neutral axis 
is enabled to shirk its duty and is largely relieved of stresses, while the 
maximum stresses will occur at points much nearer the neutral axis. On 
this theory Dr. v. Karman developed formulas by which the augmented 
deflections could be largely explained. 

Professor H. Lorenz of Danzig* offered a quite different solution based 
on the theory of thin shells and arrived at the conclusion that the stresses _ 
directly due to bending of the pipe follow the usual linear law. Deflec- 
tions calculated according to this method showed a good agreement with 
Bantlin’s experiments. 

Other investigators have used the ordinary formula for bending of 
curved bars and have made experiments which seemed to show a fair agree- 
ment with the calculated deflections. It must be noted, however, that in 
these experiments the end flanges were subject to restraint, whereby un- 
known couples were liable to be induced in the pipe bend. 

In order to determine which of the above-mentioned theories was cor- 
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rect, the author of this paper worked out a solution by a new method, as 

explained below, and applied the resulting formulas to tests on pipe bends 

carried out under his supervision at the Massachusetts Institute of Tech- 

nology. In these tests were meas- 

proressor sanrun’s ured not only the deflections of the 

srpska ends of the bends, as in former 

tests, but also the strains (stresses) 

and deformations at certain sec- 

:, tions. As in Bantlin’s tests the 

: end flanges were free to turn about 

knife edges placed normal to the 

plane of the bend. The measure- 

ments confirmed fully the theory; 

FIG. { \ the deformations, displacements 

, and stresses being in close agree- 

| ment with the calculated values. 

“A The formulas were also confirmed 

by an analysis of Bantlin’s experi- 

ments. Figure 1 shows the calculated stress distribution over one-half of 

the section of a 200-mm. pipe and for comparison the ordinary straight- 
line stress diagram is shown on the same figure. 

Consider an angular element, of angle dy, of a curved pipe bent to a 
radius R, of circular section with radius r and moment of inertia J, subject 
to a bending moment M. Let E be the modulus of elasticity. Then, 
according to the ordinary formula, the angular flexural strain ratio is given 


by 















3 2a 


Ordinal above Meutral Ais in mm 


Sthess in Lhs. per Sg. /n. 


Ady RM 
5 ae Fs (1) 


while, as shown in the Appendix, we obtain by the new method: 
Ady RM 48h?R? + 10r* (2) 
dy IE 48h?R? +r 


where h is one-half the wall thickness, 7 is the radius to the middle surface 


of the wall. 
The longitudinal direct fiber stress is found to be: 


_ Eddy 6? ; " 
P=R dy \" ~ 240R? + 5rt? (3) 


which is seen to differ from the ordinary formula for bending stresses in 
the addition of the second term inside parentheses. 
Formulas (2) and (3) were obtained after having discarded terms of 
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relatively small value and were then found to be of the same form as those of 
Dr. v. Karman. Being derived by an entirely different method and corrob- 
orated by experiments they are Pee 
believed to represent a correct ap- | | 
proximate solution. 


yf 
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y 
In calculating the deflection of a 


| 
pipe bend of the ordinary U-type | 
or of the double-offset type (Fig. 2), | 
when subject to axial forces at the H 
end flanges, we place a system of | 


rectangular codrdinates with origin | a7 





¥ S44 


FIG.2 


C2 





at the top of the bend, OX tangent 
to the axis and OY along the line - 2b, 

of symmetry of the bend. Let the 3 
entire height of the bend be H and suppose a force P is acting on each end 
flange parallel to OX tending to compress the bend. Then the bending 
moment at any point is: 
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M = PH - y) 
and the displacement of one of the end flanges, F, relative to OY is: 


aga? ys aR 48h*R? + 10r4 
Ax = Z (H — y)d(dy) = apf RU y)? “48h2R? + r* dy. (A) 


Where the bend is made up of fairly regular circular arcs, this integral 
may be evaluated algebraically, but if the bends are not true circular 
arcs it is better to find the value of the integral graphically. In the analysis 
made by the author of Professor Bantlin’s tests and of the experiments 
made at the Massachusetts Institute of Technology both methods were 
used. 


APPENDIX: DEFORMATION OF A SECTIONAL ELEMENT OF A CURVED TUBE SUBJECT TO 
BENDING 


Let R be the radius of curvature of the axis of the tube in the unstrained state, r 
the radius of the middle surface of the tube wall, and 2h the thickness of the wall. We 
assume that 2h is small relative to r and that 7 is small relative to R. The axis of the 
tube as well as the middle surface of the tube wall are regarded as inextensible. The tube 
is not subject to any pressure, external or internal. 

The tube is conceived as a body of revolution formed by revolving a transverse sec- 
tion about an axis X’X’ as indicated in figure 3. Consider a segment of the tube en- 
closed between two planes through X’X’, forming an angle dy with each other. OO’ 
is an element of the axis of the tube and PQQ’P’ is an element of the tube wall, having 
sectional end areas 2hrd@. Let p be the longitudinal stiess created by bending, then a 
force 2phrdé will be acting at each end of the element forming an angle dy with each 
other. The resultant of these forces, f, acts everywhe:e normal to the axis of revolution 
and hence normal to the neutral axis of the tube. If the bending couple is such as to 
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increase the curvature, the forces are all directed towards the neutral surface. We 
have 


f = 2hrpdody. (1) 


Due to bending, the radius of curvature R is changed to R’, dy is changed to dy + 
Ady; the ordinate y of a point P on the middle surface becomes y— Ay due to the flatten- 
ing of the section. We find the elongation 
per unit length of the filament PP’, ap- 

> proximately, 


_y Ady Ay 


Re @) 


€ 





It is important to note that f always acts 
normal to the neutral surface; it produces 
compression, shearing and bending stresses 
in the tube wall, but these stresses will not 
appreciably affect p. 

We may therefore write 


Making the elementary filament of unit 





1 
length we write dy = - 
R+y 


1 
mately R whence substituting in (1) from 


or approxi- 


FIG. 3 


(3) we get 


_2hrE(, adv _ 
omens (, dy sy) (4) 








x ee a hae! Semi elie eyo a The further development depends upon 








G . making an assumption as to the value of 
Ay. It was first tentatively assumed that Ay = ky where k is a constant over the section, 
but this led to the inconsistent result: Ay = ky*. Hence the latter value for Ay was 
adopted and led as shown below to an expression for Ay which contained also a term to 
the fifth power of y but did not differ much in value from ky’. We assume then 





Ay = ky’ = kr? sin? 6 (5) 
2hr°E/ Ad 
fu “< (ait — kr? sin? ) sin 0. (6) 


We now by usual methods calculate the deflection Ay of a transverse element of the 
tube due to the bending moment m exerted by the f-forces. By symmetry we need only 
consider one quadrant. The conditions that the sum of the angular deflections over 
the quadrant must be zero by symmetry and that the sum of the f-forces must equal 
the reaction at the neutral surface, furnish the unknown reactions. We find: 


_ bE | 


1 21 
oR _ (1—2 cos? 6) — 7 tel cost @ —6 cos?@ + a) (7) 
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The vertical deflection at any point P is found from: 


"s 
y--f 7 (e— a) ds (8) 


where i = 3h’ is the moment of inertia of the section of the tube wall. 





ie Ady . 19. ee. 
Ay = eR? [# sin? @ — kr? (2 sin’ @ + ~A sin ») | (9) 
We have 
A 
bo = (10) 
r3 


and for 6 = = we find from (9): 


rs [ady 101 22) 
Ay = — -—— 
4h*R? \ dy 120 r 


to 





which on solving for Ayg, gives: 


4 rs Ady 
Vg Oe een et 
101. dy (11) 
h2R2 sae pt 
4h?R? + 20" 


Substituting in (9) 
1 
5 cin? h2.R2 om 3 
r® sin o(s ~ ae cos °) ady 
101 dy" 
4h?R? | 4h?R? + —r' 
120 


Simplifying by the omission of small terms we may write: 





Ay = 


6r' sin? @ = =©Ady 
Ay a ee a 12 
y * 2408R? + br dy (12) 


which is consistent with the assumption made in (5). Substituting in (3) from (12) 


we find the stress: 


me. REE cr AM Ga 
Le E(s ane) dy (18) 


In order to determine the angular deflection we equate the external bending moment 


to the moment about the neutral axis of the stress forces: 


um [amp o- aye 
0 


Ad ; 
and neglecting terms in (“ ) we find on solving for the angular deflection: 


Ady _ MR 48h*R? + 10rt (4) 
dy El 48R?+ fr 


We are now able to find the stress at any point from (13) and the displacement of the 


ends of a pipe bend from formula (4) in the main body of the paper. 
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1 A more detailed account will appear in J. Math. Physic., Mass. Inst. Tech., Vol. 6, 
No. 1 (1926). 

2 Mittheilungen tiber Forschungsarbeiten, Heft 96, 1910. 

3 Dingler’s Polyt. Journal, 1912. 

4S. Crocker and S. S. Sanford, Amer. Soc. Mech. Eng., 1922. 


THE HEAT CAPACITY OF SOLID ALIPHATIC CRYSTALS. II 
By E. O. SaLant* 
Jouns Hopkins UNIVERSITY 


Communicated May 10, 1926 


In the preceding issue of these PROCEEDINGS, the fact of characteristic 
frequencies for non-polar bonds was introduced into the Born-Karman 
heat capacity equations, which were then, with certain approximations, 
brought into a general simplified form for non-polar compounds.'! With 
the aid of these simplified equations, the heat capacities of some organic 
solids will now be calculated and compared with their observed heat 
capacities. 

Due to the still scanty knowledge of the infra-red vibrations of organic 
substances, it will be necessary to make use of approximations to get the 
required frequencies, and so agreements and disagreements of the calcula- 
tions with the measurements cannot be taken as severe tests of the validity 
of the equations. 

For the molecular heat capacity at constant volume, C,, of an open- 
chain compound A,B,C, having m atoms and } non-polar bonds to the chem- 
ical molecule, we had 


C, = f'p (m/T) + 2[(x—-WE’ (v/T) + yE’ o/T) + 2E’ (vX/T)] + 
DEW /T) (1) 


Im 


where f'n (v,,/T) is the Debye atomic heat capacity equation and E’(v/T) 
is the Einstein function for a linearly oscillating atom. These are defined 


as 
__9R_ %(0/T)'¢"7 de | 
f'v Um/T) = won [rae | (2) 
: = R(6/T)*e’". » hv 
e OY aie Fs 


where JT is the temperature, h is Planck’s constant, R the absolute gas 
constant and k the constant per molecule (R divided by the Avogadro 
number). 
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Equation (1) will now be limited to the case of aliphatic compounds of 
C, H and O. It is necessary to choose first some atom as nearest to the 
center of the molecule, that is, some atomic species corresponding to the 
species A. For the substances considered in this paper, a C atom will be 
chosen; this may not always be proper (for example, if methyl or ethyl 
ether were to be considered, the O atom would have to be the central atom). 
We have, then, for the compound C,H,0, 


Cy = f'o(vm/T) + U(x — IE’ (¥©/T) + yE! 0 /T) + 
b 
2k! (v /T)] + DUE" /T) (3) 


where »,, refers to the molecular vibration »©, v™ and » to the vibra- 
tions of the C, H and O atoms across the bonds, and the last term to the 
bond vibrations. Equations (2) and (3) permit the calculation of the heat 
capacities of aliphatic solids if the required frequencies are known; methods 
of obtaining values for these frequencies will now be discussed. 

In the first attempt to apply the quantum equations of heat capacities 
of solids to aliphatic compounds, the writer found it necessary to use a 
mean vibration in all directions for the atoms and to obtain this frequency 
with the Koref-Lindemann melting-point formula.? If the atoms are 
thought of, first as vibrating with the frequencies given by this formula, 
but then restricted to the characteristic bond frequencies along the lines 
of the bonds, and the effect of the bond on the other vibrations is neglected, 
then this formula may be used for calculating the vibrations of the atoms 
across the bonds. Since there is no method as yet for obtaining those 
frequencies (infra-red spectra being still inadequate), the Koref-Lindemann 
formula had to be used in this way for the 2(x — 1 + y + 2) frequencies. 
The writer regards this as the greatest source of error in these calculations. 

The mean frequency v,, for the Debye equation is determined in the 
strict theory by the elastic properties of the crystal, but measurements of 
elastic constants of organic compounds are, unfortunately, lacking. For 
those compounds whose heat capacities have been measured at the lowest 
temperatures (glucose, for example, at 20°K.), where the Einstein functions | 
become negligible, v,, can be obtained from the heat capacity measured 
at some one of these lowest temperatures. Again, unfortunately, too few 
organic heat capacity measurements have been carried to these low tem- 
peratures. 

The remaining alternative is to use the melting-point equation, following 
Nernst,’ with the molecular weight M in place of the usual atomic weight: 


Ym = 3.1 X 10"7T,d" un’ 


where Ts is the absolute melting-point and d is the density of the com- 
pound. 
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For our purposes, this may be approximated still further. ‘The density 
of organic compounds being usually not far from 1, we may write 


Vm = 3.1 X 10°T,“M-* 


which enables us to evaluate the Debye equation for all those organic 
compounds whose molecular weight and melting-point are known. For 
comparison, we have the molecular heat at a very low temperature of glu- 
cose only: at T = 20, from the melting-point formula and Debye equation, 
c = 0.18, measured 0.09. 

The characteristic frequencies of the C—H, C—C, C—O and O—H 
bonds are required for our calculations. The small variations of these 
frequencies with variations in the degree of saturation and loading of the 
atoms of a bond will be neglected. 

The absorption band of the C—H bond was found at 6.94, of the C—C 
at 28.0y, both by Ellis.‘ 

The absorption of the alcoholic OH group was found at 2.954 by Cob- 
lentz,> at 3.0u by Weniger.* Coblentz further demonstrated that this 
band disappears in organic acids, and neither in his nor in Weniger’s 
measurements of the acids can any band characteristic of carboxylic OH 
be detected. ‘The band at 3y will, therefore, be taken as the fundamental 
for the alcoholic O—H bond. 

Coblentz’ lists a band at 29.4u as characteristic of carbonates. The 
measurements of Coblentz and of Schaefer and Shubert*® showed a band 
varying from 13.32 to 15.204 for various carbonates, and this may be 
taken as the approximate second harmonic of the 29u band as a funda- 
mental (first harmonic). 

If we assume that these are due to the C—O bond, evidence should ap- 
pear in the spectrum of CO,. Rubens and Aschkinass® found an emission 
band of CO, at 14.1, and absorption at 14.74; Burmeister!® showed the 
absorption to be a double band, 14.70 and 15.05n. 

Weniger attributed an absorption band at 5.94 to the C=O group, 
and Marsh" gives further evidence of a C=O vibration at 5.8u. It is 
to be noted that this would be a fifth harmonic of a fundamental 294 C—O 
band. 

Weniger found a band from 9.1 to 9.84, common to alcohols, esters and 
organic acids, which he could not classify (this being before the assignment 
of bands to bonds, initiated by Ellis). The one bond common to all of 
these (in addition to the C—H with band at 6.9) is the C—O, and a band 
in this region is, therefore, to be looked for in ethers. Coblentz’s curves!” 
show absorption bands at about 9.34 for both methyl and ethyl ether 
vapors. The band at 9.1 to 9.8u characteristic of the C—O bond would be 
an approximate third harmonic of the fundamental 29y. 

Thus this fundamental serves to correlate bands in CO, carbonates, 
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ketones, ethers, alcohols, acids and esters, and will, therefore, be taken as 
the characteristic of the C—O bond. 
The hv/k,values of these bonds are: 


CE c-6 co O—H 


2115 521 503 4863 


It is to be noted that the contributions to the heat capacity of the 
O—H bond are inappreciable below 347°K. For the carboxyl O—H, 
in formic and oxalic acids, the H atoms are given a mean vibration in all 
directions, calculated by the Koref-Lindemann formula. 

Thus the mean atomic heat at constant volume, c,(—C,/n), was ob- 
tained by equations (2) and (3), the melting-points, molecular weights and 


GLYCEROL Oxa.ic AcID 
Cp Cp 
¥ CALC. ops. 4 T CALC, oss. !8 
280 3.09 2.46 ; 200 3.18 2.30 
227 2.54 1.79 160 2.63 1.95 
174 1.97 1.44 120 2.00 1.60 
116 1.25 1.07 100 1.60 1.41 
70 0.69 0.72 
GLUCOSE Formic AcIp 
Cp Cp 
x CALC. oss,'7 ¥ CALC. oss.!5 
300 3.06 2.29 270 3.60 3.42 
260 2.69 1.95 237 3.39 3.02 
220 2.34 1.64 5 205 3.18 2.73 
180 1.93 1.34 176 2.97 2.48 
140 1.40 1.06 94 2.06 1.73 
100 0.84 0.78 71 1.57 1.45 
60 0.34 0.49 
20 0.18 0.09 
MeEtuyt ALCOHOL Etsy, ALCOHOL 
Cp Cp 
T CALC. oss.'3 : CALC, oss. '5 
152 2.09 2.36 140 1.78 1.81 
112 1.64 1.91 130 1.64 1.60 
101 1.53 172 120 1.51 1.51 
89 1.27 1.61 110 1.39 1.42 
100 1.26 1.31 
90 1.15 1.25 
Butyt ALCOHOL ACETONE 
cp Cp 
T CALC. oss. 18 sy CALC. oss. '!6 
151 1.63 1.35 173 2.32 2.32 
145 1.52 1.32 153 2.05 2.09 
111 1.10 Be 133 1.74 1.80 
100 0.96 1.03 113 1.45 1.64 
91 0.85 0.95 93 1.17 1.22 
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characteristic bond frequencies. The mean atomic heat at constant 
pressure, c,, was then calculated by the Nernst-Lindemann relation:’® 


Cy — & = oT A,/Ts 


where A, = 0.0214. This, again, had to be used in place of the thermo- 
dynamic relation, due to the lack of measurements of thermal expansions 
and compressibilities of organic solids. 

Tables of the calculated and observed heat capacities of methyl, ethyl 
and butyl alcohols, glycerol, acetone, glucose and formic and oxalic acids 
are given, from the lowest to the highest temperatures measured. 

Comparison of the results of this paper with those of the preceding one on 
aliphatic heat capacities? shows that the calculated values of the three 
alcohols, of acetone and of formic acid are now uniformly nearer the 
experimental values, the improvement at lower temperatures being es- 
pecially marked. At these temperatures, the calculated heat capacities 
of glycerol, glucose and oxalic acid are also better, but at higher tempera- 
tures are now worse. However, in the calculations of that paper a special 
assumption as to the C vibrations in these three substances was made; 
without this special assumption, their values even at higher temperatures 
are now improved (for example, error in glucose at 300°K. formerly 
50%, now 30). 

It may be concluded, then, that by introducing asymmetric vibrations 
for the atoms due-to their bonds, and also molecular vibrations, the calcu- 
lated heat capacities of these organic compounds show consistently better 
agreement with measurement. These and other results are to be ‘more 
fully discussed at a later date. 
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THE EFFECT OF CRYSTAL SIZE UPON THE INTENSITY OF 
X-RAY REFLECTION 


By R. J. Havicnurst! 
CruFt LABORATORY, HARVARD UNIVERSITY 


Communicated May 3, 1926 


It appears possible to determine with considerable accuracy the electron 
distribution in atoms? from a knowledge of the intensity of reflection of 
X-rays by crystals containing those atoms. However, the effect of cer- 
tain other factors, notably the size of the reflecting crystal, upon the in- 
tensity of reflection must be known before conclusions can be drawn 
concerning electron distribution. 

Darwin’ has studied theoretically the effect of crystal size upon the in- 
tensity of X-ray reflection. His expression for the intensity of reflection, 
I, at an angle 6 from the face of a small perfect crystal may be written 


(1 + cos? 26) _»*sin?0 (1 + cos? 26) 

panies a wi sin20. Q. (1) 

N is the number of atoms per unit volume; f is a function of _— and 
of the electron distribution in the atoms; os is the Debye tempera- 
ture factor. If the crystal has a linear dimension larger than 10~° cm., 
the intensity of reflection is modified by ‘‘extinction.”” Before considering — 
extinction, let us define a perfect crystal as one in which the planes of each 
unit cell are perfectly parallel to the corresponding planes of every other 
unit cell, so that the whole crystal can be regarded as a macroscopic copy 
of the unit cell. Darwin distinguishes two kinds of extinction, primary 
extinction, which occurs within each little block that acts like a perfect 
crystal, and secondary extinction, which represents the shielding of the 
blocks deep within the ‘crystal by the reflection of the X-rays from the 
blocks near the surface. The effect of primary extinction is represented 
by the correction factor in the expression 


tanh 2Q . ~ “ig 
0. = 0 — ee (2) 
\ d? cot 9 











20-> 


Q’ is proportional to the measured intensity of reflection and becomes 
equal to Q in (1) when extinction is negligible; d is the thickness of the 
crystal. If expression (2) is evaluated for rock-salt, using the value of 
Q’, obtained by Bragg, James and Bosanquet for the (200) reflection‘ 
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and corrected for secondary extinction, as the Q under the radical, it is 
found that ford = 5 X 10-*cm., the correction factor is 0.96, while for d = 
5 X 10-* cm., the factor is 0.03. Since the value of Q used to obtain these 
figures is really a Q’ value and may have been reduced by primary extinc- 
tion, it is evident that these are maximum values for the correction factor. 

Secondary extinction is appreciable only in crystals larger than 10-?- 
10-* cm. and is directly proportional to the intensity of reflection. Bragg, 
James and Bosanquet have been able to correct satisfactorily for its effect 
in their measurements upon large single crystals of rock-salt. Darwin 
suggested that information regarding extinction might be obtained by meas- 
uring the intensity of reflection from crystal powders of varying sizes. 
If each powder particle is a perfect crystal, the intensity of reflection would 
vary with the size until the particles were less than 5 X 10~° cm. in thick- 
ness. Further, the relative intensities for different planes would vary, 
because the extinction depends upon the intensity and angle of reflection. 

Evidence which has been gathering from different sources® indicates that 
in spite of its flawless appearance, a macroscopic crystal is almost never 
perfect, but is a mosaic structure consisting of many. very small perfect 
crystals with very nearly but not exactly the same orientation. The size 
of these perfect crystal blocks and their divergence from the mean orienta- 
tion depend upon the composition of the crystal and upon the mechanical 
and thermal treatment which it has undergone. Mark estimates the 
linear dimensions of the perfect crystal blocks in most crystals to be about 
10-4cm. Smekal estimates the thickness of the perfect crystals in rock- 
salt to be 4 X 10->cm. Becker obtains a value of less than 10—° cm. for 
tungsten. 

That there is some primary extinction operating to reduce the intensity 
of reflection from a large single crystal is evidenced by the fact that grind- 
ing a crystal face increases the intensity of reflection to a marked degree. 
Mechanical strains probably break up the crystal blocks near the surface 
or produce distortions in them which reduce primary extinction. Grind- 
ing a crystal to a powder, then, might be expected to distort the powder 
particles so that the size of the perfect. crystals would be many times 
smaller than the size of the particles themselves. The experiments de- 
scribed below were carried out in an attempt to decide whether the re- 
flection of X-rays by powdered crystals is free from the effects of extinction. 

Experimental—The relative intensities of reflection of the Mo Ka 
doublet by the different planes of many crystal powders were determined, 
using the radiation, filtered through enough ZrO, to reduce Mo K(7 + 8) 
to 0.4% of Mo K(a: + a), from a Mo tube operated at a constant po- 
tential of 35,000 volts. The primary beam was reflected from the face 
of a briquet of the powder under investigation into an ionization chamber, 
the distances from focal spot to powder, a, and from powder to ionization 
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sina 4a 
sin B b 
tween the primary beam and the surface of the powder and 8B = 20—a. 
This arrangement, which has been described fully by Brentano,’ causes 
a sharp focussing of the reflected beam, even though a primary beam of wide 
divergence is used. With this apparatus, it is possible to measure the in- 
tensities of reflection out to@ = 45°. Intensity measurements on powders 
may be compared with those on single crystals by the application to the 
I sin? 6 cos 6 
j 1 + cos? 26 
measured intensity of reflection and 7 is the number of planes belonging to 
the form which produces the reflection. Absolute measurements of the 
reflecting power have not been made. 

Sodium Chloride.—Bragg, James and Bosanquet have measured the 
intensity of reflection of Rh Ka rays from single crystals of rock-salt, 
and have corrected their re- 


chamber slit, }, determining the ratio 





where a is the angle be- 


powder method data of the expression A? = where I is the 





sults for secondary extinc- NaCl 
tion. Their method was un- 
able to take account of pri- 

. . me NES 4-BJBlUncorrected) 
mary extinction. Theupper “TF . PAL ot em 
and lower curves in figure 1 o-deyte ot 
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plotted against (h?+ k?+/*), “| 
The relative values of A? for 4 
the different planes are in- 
dependent of wave-length, so 
the powder measurements 
made with Mo Ka radiation 
may be plotted in the same 
figure. Since primary extinc- 
tion is proportional to the _,, 
intensity and to the angle 
of reflection, any variation 
in extinction should be evi- i, 
denced by a change in the 
relative values of A*. Sec- 
ondary extinction, which is proportional to the intensity of reflection only, 
should also be detected by a change in relative intensities, as is shown 
by the difference between the corrected and uncorrected curves of Bragg, 
James and Bosanquet. 

Sample 3 was a powder which had been ground until the particles were 
all smaller than 4 X 10-*cm. The intensity measurements are seen to fall 
upon a smooth curve which coincides with that representing the data of 
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Bragg, James and Bosanquet, when the latter are corrected for secondary 
extinction. Some of the same powder was ground mechanically for 
several days, until the particle size was 1 X 10-*—2 X 10% cm., and gave 
results identical with those of 

Ca F; sample 3. Sample 6, a powder 
‘ which had passed a 200-mesh 
sieve and been held back by 
a 325-mesh sieve, with particle 
size 4.4—7.4 X 107% cm., was 
heated at 760° for 40 hours, 
then cooled gradually to room 
temperature. Because of the 
high vapor pressure at 760°, 
the crystals had grown. After 
having been reground to the 
original size, they gave results 
in agreement with sample 3. 
Samples 7 and 8 were from the 
same lot as 6 but were heated 
at a lower tethperature, 700°, 
for 16 and 46 hours, respec- 
tively, then cooled slowly to 
room temperature. The part 
of the powder which remained 
Se ~«C Within the original limits of 
(h*+ k* 1? size was investigated, care 
wt Saat being taken to avoid mechan- 
ical strains. The results for these two samples were in close agreement, 
and are plotted in the figure. Evidently sufficient diffusion occurred within 
each powder particle to relieve some of the distortions, and a certain amount 
of primary extinction took place. Analogous results have been obtained by 
James® onheating single crystalsof rock-salt. The powder particles were too 
small to permit of any noticeable effect from secondary extinction. Sample 
8, after being ground for two hours, was reinvestigated, and then gave re- 
sults in agreement with those from sample 3. Sample 9 was prepared by 
condensing the vapor of NaCl which had been vaporized from a hot plati- 
num wire. Most of the particles were mere points when viewed through 
the oil immersion lens of a microscope, about 10~* cm. in diameter, but a 
few were as large as 3 X 10-*cm. ‘The results show a slight amount of 
primary extinction, which modified the most intense reflection. Probably 
crystals formed from the vapor phase are more nearly perfect than those 
formed from solution, and even though most of the particles are so small 
that primary extinction is negligible, the few larger particles produce an 
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extinction the effect of which is measurable. Although sample 9 was 
later ground for two hours, there was no appreciable change in the inten- 
sities of reflection, possibly because the particles were so small as to undergo 
no strain during the grinding. 

Calcium Fluoride——MacInnes and Shedlovsky® have measured the rela- 
tive intensities of reflection of Pd K rays from various planes of fluorite, 
but have made no attempt to correct their data for secondary extinction. 
Their data is plotted along with the powder method measurements in fig- 
ure 2. The upper curve represents reflecting planes in which fluorine re- 
enforces calcium, while the lower curve represents reflections from cal- 
cium planes alone. It is evident that extinction has here produced a 
notable effect. Three different samples of powdered CaF», all giving the 
same results, were used. Sample 1 was made by grinding a clear colorless 
fluorite crystal to 4.4—7.4 X 10-*cm. Sample 2 was from the same crys- 
tal but was further ground for 1 day to 4 X 10-*—1 K 10-‘cm. Sample 
3 was a pure precipitated powder with particles all smaller than 4 x 107% 
cm. It is possible that in addition to secondary extinction, which most 
certainly affected the results of MacInnes and Shedlovsky, primary 
extinction was also effective. As evidence to the contrary, however, it 
should be pointed out that both the (111) and the (400) reflections, which 
MacInnes and Shedlovsky found to be equal, are equal in the powder 
method results. Since primary extinction depends upon the angle of re- 
flection as well as upon the intensity, the correction for reflections at the 
two different angles represented by (111) and (400) reflections should only 
be equal in case secondary extinction alone was operative. 

Calcite —Although it is difficult to make powder method measurements 
because of the great number of reflecting planes, measurements have 
been made with sufficient accuracy upon samples of ground calcite and 
precipitated calcium carbonate to show that the more intense reflections 
measured by W. H. Bragg’ are several hundred per cent in error because of 
extinction. 

Conclusions.—The conclusion to be drawn from this work, along with 
that of other investigators, is that primary extinction is absent in crystals 
that have been grown from solution and ground to less than 10~* cm. in 
thickness and that measurements upon the intensity of reflection of X- 
rays by these powders may be used in the determination of electron dis- 
tributions. Further, in crystals as poor as the alkali halides, single crystal 
measurements similar to those of Bragg, James and Bosanquet are satis- 
factory after correction for secondary extinction. 
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THE INTENSITY OF REFLECTION OF X-RAYS BY LITHIUM, 
, SODIUM AND CALCIUM FLUORIDES 


By R. J. Havicuurst! 
CruFtT LABORATORY, HARVARD UNIVERSITY 


Communicated May 5, 1926 


The scattering power for X-rays of an atom in a crystal falls off with 
increasing angle of scattering in a manner that is determined by the elec- 
tron distribution in the atom. In an effort to determine the electron dis- 
tribution in certain light atoms, X-ray measurements have been made on 
the scattering power of powdered crystals of LiF, NaF and CaF». 

Darwin's expression? for the power P diffracted in a cone of semi-apex 
angle 26 by a bit of crystal powder of volume 5V may be written 

i a as ite OO a ae 
I is the power in the primary beam; 2 is the number of atoms per unit 
volume; j is the number of planes in the form under investigation; F? 


is a function of — and of the arrangement of the electrons in the atoms, 
2 sin? 

and includes the Debye temperature factor r. F, which is called 
the ‘‘atomic structure factor,’’ approaches the number of electrons in the 
atom at small angles of reflection and theoretically should become equal 
to this number when sin@ = 0. All of the unknown quantities in the above 
expression except F may be measured, so that F may be calculated. Since 
the measurements reported in this note are only relative, we may dispense 
with the part of (1) which does not enter directly into the determination of 
the relative intensities. We then have 
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Paine 
ae 7 1+ cos? 26 (2) 
The power is measured by an ionization chamber at a distance b from the 
powder sample, with a slit of length / which is small compared with b 
sin 28. If the measurement is made so as to take in all of the angles in the 
neighborhood of 2@ at which any measurable power is diffracted, the 

, tere ‘ Pl 

power entering the ionization chamber is P, = onb an 26 
Expression (2) then becomes 


P, sin? 6 cos 6 


Big 4? Se. 
44 7 1+ cos? 20 


(3) 


In order to get sufficient power in the scattered rays for accurate meas- 
urement at large angles, a focussing method described by Brentano* has 
been used. The primary beam from a Mo tube, filtered through enough 
ZrO, to reduce Mo K(y + 8) to less than 0.4% of Mo K(a; + as), was 
reflected from the face of a briquet of the powder under investigation, the 
distance from focal spot to powder, a, and from powder to ionization cham- 
a 
b . 
primary beam and the surface of the briquet, and 8 = 20—a. The tube 
was operated by a high potential storage battery‘ at 35,000 volts. 

Relative intensities could be determined to an accuracy of considerably 
less than five per cent except in the case of the very weak reflections. 
The conceivable sources of systematic error have been investigated. It 
was found that the pressure used in making the briquet has a tendency to 
orient the surface layer of crystals if they are larger than 10~* cm. in thick- 
ness. Consequently, the original surface of each sample was shaved off. 
The effect of particle size upon the intensity of reflection has been in- 
vestigated (cf. preceding note) and it is believed that extinction was not 
operative in these experiments. 

In figures 1 and 2, F curves for the different classes of crystal planes and 
their component atoms are plotted against sin @. Since we are here dealing 
with relative values of the atomic structure factor, the ordinates will be 
called F’. However, an effort has been made to place these F’ values on 
an approximate absolute scale. The author’s F’ curve for Na in NaCl 
was evaluated with the aid of the data of Bragg, James and Bosanquet 
on rock-salt,> and an F curve for Na was thus obtained and compared with 
the F’ curve of Na from sodium fluoride. These two curves for Na could 
not be made to coincide, presumably because the electron distribution in 
the Na atom of NaCl is different from that of the Na atom in NaF. The 
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values of F for sin @ = 0.20 on the two curves were arbitrarily made 
equal. F’ values for NaF were then determined and since all the crystals 
have fluorine in common, the F’ curves for fluorine were made to agree 
as closely as possible. A two-fold assumption is therefore the basis of the 
valuation of these F’ curves: 
(1) the scattering power of Na 
in NaCl and NaF is the same 
when sin 6 = 0.20; (2) the 
scattering power of fluorine 
in LiF, NaF and CaF, is ap- 
proximately the same at all 
angles. We may now consider 
the curves of figure 1. 
Lithium Fluoride—The 
sample was carefully freed 
from Ca and Mg and several 
times reprecipitated. The 
upper curve represents planes 
with all even indices, in which 
Li and F reénforce each other; 
in the lower curve, for planes 
with all odd indices, Li inter- 
feres with F. Debye and 
Scherrer® have determined the 
relative intensities of reflection 
of the first nine planes of LiF 
by photometric study of pow- 
der method photographs. 
“sing®* °F 96 — °1 Their results fall reasonably 
eS well upon the curves of figure 1. 
They interpreted their work as proof that the crystal lattice consists of Lit 
and F~ ions, rather than of neutral atoms. We have seen that the F values 
should be equal to the number of electrons in the atoms composing the 
crystal planes at zero angle of reflection. Consequently the F + Li curve 
should approach the value 12, and the F — Li curve should approach the 
value 8 if the lattice consists of ions, or 6 if the lattice consists of neutral 
atoms. The corresponding ratios of ; - 7 are 1.5 and 2. Debye and 
Scherrer, by extrapolating the curve representing this ratio to zero angle of 
reflection, obtained a value of approximately 1.5. While there is abundant 
evidence from other sources that the lattice of the alkali halides is ionic, 
the shape of the F curves at small angles is entirely unknown to us and re- 
liance can hardly be placed upon curves extrapolated into this region. The 
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. F+ Li 
ratio Foi : 
ing the smaller value at small angles. 

Sodium Fluoride —Samples from two different sources, each sample con- 
taining about 99% NaF, were investigated, with results which were in satis- 
factory agreement. The curve labelled Na—F is of particular interest be- 
cause it represents the difference in scattering power of the Na* and F- 
ions, on the assumption that the lattice is ionic. This curve should pass 
through zero at sin 6 = 0. Since the ions have the same number of 
electrons, there seems at first to be no reason for assigning to Na* the larger 
scattering power. But Nat, possessing a nuclear charge of +11, would 
bind its electrons more tightly than F~, with a nuclear charge of +9. 
A system with electrons nearer its center has a larger scattering power 
through the angular domain covered by the curves of figure 1. 

Calcium Fluoride—Samples of ground fluorite, and a precipitated prod- 
uct containing less than half of one per cent impurity, gave similar results. 
Fluorite, with a crystal structure different from that of the alkali halides, 
gives reflections which fall into three classes, as shown in the figure. The 
Ca—2F reflections are very weak, and occur so close to the Ca reflections 
that their measurement is quite difficult. Only two direct measurements 
were made, shown by the circles on the curve; they are correct only in order 
of magnitude. The points marked + on the two lower curves were 
obtained by an indirect method. The combined intensity of a Ca reflec- 
tion and its neighboring Ca—2F reflection was measured. ‘Then an inter- 
polated value from the Ca curve was used for the intensity of the Ca re- 
flection and the amount left over from the total was assigned to the Ca—2F 
reflection. ‘These values also are only approximate. The results of Mac- 
Innes and Shedlovsky’ for the relative intensities of reflection from fluorite 
are affected by extinction (cf. preceding note), but the magnitudes of the 
weaker reflections have been located approximately in the figure, and 
marked X. ‘The broken curve is that obtained for Ca—2F from the two 
upper curves. Ifthe Catt+ ion, with 18 electrons, is assumed to be present 
in the lattice along with 2F~- ions, with 10 electrons apiece, the curve 
representing Ca—2F must drop to zero and then rise as a 2F —Ca curve toa 
value of 2 when sin 6 = 0. Here we see that Ca**, with an excess posi- 
tive charge, has drawn its electrons in and consequently scatters X-rays 
more strongly than the two F~ ions, which possess an excess of negative 
charge. 

Atomic Scattering Curves.—In figure 2 are the atomic scattering curves 
obtained from the curves of figure 1. Of considerable interest is the fact 
that the curves for fluorine from different compounds are nearly alike. 
Provided that these curves do not differ on an absolute scale, the conclu- 
sion is that the fluorine ion is in a force field of the same magnitude in all 


in figure 1 varies between 2 and 1.5, and seems to be approach- 
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three crystals. As has already been pointed out, the Na curves from NaF 
and NaC! do not coincide, probably because the chemical forces are of dif- 
ferent magnitudes in the two crystals. Examination of the heats of forma- 
tion of the compounds in question gives qualitative support to this ex- 
planation. ‘The fluorine curve of CaF, possesses a slight hump, as do the 
Ca and Licurves. If experimental error is ruled out, the humps indicate 
either a definite shell-like structure of the electrons, or a lack of spherical 
symmetry in the atoms. Lit, with only two K electrons, could: hardly 
give anything except a smooth scattering curve unless the system lacks 
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FIGURE 2 


sphericalsymmetry. The electron distribution curves which are now being 
worked out for these crystals by a method involving a Fourier analysis® 
should help to answer many of the questions raised by the appearance of 
these curves. 

The atomic scattering curves show clearly that the assignment of scat- 
tering power to atoms in accordance with their total number of electrons 
is only a crude approximation. A similar variation of scattering power 
with angle of scattering is predicated for all ions by such a procedure. 
Especially uncertain is this assignment of scattering power in the case 
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of oppositely charged ions (the usual case in crystal analysis) because the 
excess of positive or negative charge produces a tightening or loosening of 
the electron atmosphere of the original atom which profoundly affects its 
scattering properties. 
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NOTE ON THE POSTULATES OF THE MATRIX QUANTUM 
DYNAMICS 


By J. H. VAN VLECK 


DEPARTMENT OF Puysics, UNIVERSITY OF MINNESOTA 
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Communicated April 22, 1926 


Through the researches of Born, Heisenberg and Jordan,'** and of 
Dirac,**® the dynamics of the quantum theory have been formulated in 
terms of matrices. The fundamental postulates made by Born, Heisen- 
berg and Jordan are the Ritz combination principle 


re Se a ee ee 


v(nm) + v(mk) = v(nk) ; (1) 
the Hamiltonian canonical equations® ' 
an : | 
Gq: = Op, Py hr oq, (2) 
the quantum conditions, } 
h i 
Ped — GPx = 51 (3) 
and the commutability relations 
- = 0, 1#k), 

Pde — OePr ( ) (4) 


de — 4G: = 9, PP — PrP: = O. 


Here H is the energy (Hamiltonian function) and the q’s and p’s are 
coordinates and momenta. We shall suppose that there are s degrees of 
freedom, so that the subscripts k and / range from 1 tos. Allexpressions 
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printed in bold-face type are matrices, as we follow Born’s notation through- 
out. From (1), (2), (3) and (4) Born, Dirac and others deduce the Bohr 
frequency condition 


hv(nm) = H(nn) — H(mm) (5) 


and the conservation of energy _ 
H = 0. (6) 


The purpose of the present note is to call attention to the fact that (2), 
(4), (5) and (6) can be taken as the basic postulates of the new theory 
instead of (1), (2), (3) and (4). The fundamental postulates are then the 
canonical equations, the commutability of the order of multiplication in all 
the products q,q), PxP;, GzP, except q,p,, the Bohr frequency condition, and 
the conservation of energy. In other words, if we accept (2) and (4), the 
equations (1) and (3) are not merely sufficient (as proved by Born, Heisen- 
berg and Jordan*) but also necessary for the validity of the Bohr frequency 
condition and the conservation of energy. This is not surprising, but is per- 
haps worth noting inasmuch as it does not appear to have been specifically 
mentioned in the various papers previously written on the matrix dynamics. 
In fact, the Bohr frequency condition and the conservation of energy are 
more familiar concepts than the quantum conditions (3) and so may ap- 
peal to some readers as more lucid postulates than (1) and (3). One must 
still make the postulate (4) that the order of multiplication is commuta- 
tive in all products of the form 4q,q), p,p;, Q,P), except q,p,, and un- 
fortunately the basis for this assumption is not apparent unless, following 
Dirac,* we establish a correspondence with the Poisson’s bracket expres- 
sions of ordinary dynamical contact transformation theory. In systems 
with one degree of freedom, however, there is only one q and p, and one can 
then dispense with (4) entirely, so that the fundamental postulates take 
the very simple form (2), (5) and (6). 

Proof —The Ritz combination principle (1) is an obvious consequence 
of the Bohr frequency condition (5). To show that postulates (1) and (3) 
can be replaced by (5) and (6) we must in addition deduce (3) from (2), 
(4), (5) and (6). Our proof will be very similar to Born, Heisenberg 
and Jordan’s converse derivation® of (5) and (6) from (1), (2), (3) and (4). 
By virtue of the Bohr frequency condition and the constancy of energy, 
we have § 

ha, = 2xi(Hq, — 9H). (7) 


This is a well-known result. It follows directly inasmuch as an element 
of the matrix on the left-hand side of (7) is 


h.2riv(nm)a(nm)exp.2riv(nm)t 
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while an element of the right-hand side is 
2ni(H(nn) — H(mm))a(nm)exp.2riv(nm)t 


in consequence of the rules for multiplication of matrices and the fact that 
H is a diagonal matrix because of (6). Here a(mm)exp.27riv(nm)t denotes 
an element of the matrix q,. Equation (7) is thus equivalent to (5) 
whenever the q, and p, can be represented by matrices, as in multiply 
periodic systems. Equation (7) is, however, more general than (5) inasmuch 
as (7) is applied by Born and Wiener’ and by Dirac’ even to “‘aperiodic’”’ 
motions which cannot be represented by ordinary matrices. The derivation 
which we shall give of (3) applies to aperiodic systems if we take (7) as a 
fundamental postulate in place of (5) and (6). 
Comparing (2) and (7) we see that 


H H = = S 8 
Aa oe (8) 
We shall next establish the equation 
of h 
i ES Pe ae (9) 


where f is any function of the p’s and q’s. This equation was deduced by 
Born, Heisenberg and Jordan from (1), (2), (3) and (4), whereas we must 
deduce it from (2), (4) and (7). Before proceeding to the proof we must 
define what is meant by the derivative of one matrix with respect to 
another, as such a derivative is involved on the right-hand side of (9). 
We shall assume that the operation of matrix differentiation is defined in 
such a way that 


Op: _ oe ae 

a a ee ss (10) 
ra) Oa ob re) Oa ob 
am (t. +b). oe bees Oe Rhee 
Set eae a ae ee 


where a and b are any two functions of the coérdinates and momenta. 
This procedure is tantamount to taking the derivative to be a ‘‘differential 
quotient of the first kind’ in the terminology of Born, Heisenberg and 
Jordan.* In fact the “differential quotients of the first kind’”’ are perhaps 
more readily understood if defined by (10) rather than by the alternative 
limiting process used by Born, Heisenberg and Jordan.* 
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Now by (4) and (8), equation (9) is obviously valid in the particular case 
that f is any one of the 2s expressions 


Pi,---» De—1» Petis---» Dsy qi,.--, q;, H. (11) 


Consequently (9) is valid when f is any function of the 2s arguments (11). 
To establish this we need only note, with Born, Heisenberg and Jordan,* 
that functions of matrices are defined by successive applications of the 
operations of addition and multiplication, and that if (9) holds forf = a 
and f = b, it alsoholdsforf = a+ bandf=ab. Inthecaseoff=a+b 
the truth of this statement is obvious, while for f = ab we have (cf. 
Born’) - 


abq, — q,ab = a(bq, — q,b) + (aq, — q,a)b 


db da h r) h 
wy BOS... NG See ae — by (10). 
(a sigs =» |= ~ (ab) 5 y (10) 


Since (9) thus holds for any function of the 2s arguments (11) it also 
holds for any function of the 2s coérdinates q; and p,; (J = 1,...s). This 
follows since by solving the energy equation 


H = H(q,. - 9s, Di,-- -Ps) 


for the variable p,, we can express p, in terms of the arguments (11). 
Thus p;, can be eliminated and the arguments of f always regarded as (11). 

As we have shown (9) holds for any function f, we can, in particular, 
take f = p,. This gives immediately the desired result (3), Q.E.D. We 
could also easily establish the relations 


of h 
itty ee te 


which are analogous to (9) but which are not needed in proving (3). 
The writer wishes to express his thanks to Prof. Tate for suggesting the 
question discussed in the present note. 
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CONFORMAL INVARIANTS 
By JosePH MILLER THomas* 
PRINCETON UNIVERSITY 
Communicated April 1, 1926 


1. Ihave previously! pointed out the existence of a conformal invariant, 
the conformal connection, which in the problem of the conformal corre- 
spondence of Riemann spaces plays a réle similar to that of the Christoffel 
symbols in the problem of the equivalence of two quadratic differential 
forms, and by expressing integrability conditions of its law of transformation 
given a derivation of the conformal curvature tensor discovered by Weyl. 
By continuing the process of elimination employed in the former paper, I 
now show how a complete set of integrability conditions for the equations 
of transformation of the conformal connection can be obtained. These 
integrability conditions express the laws of transformation of a set of con- 
formal invariants. In addition to the first partial derivatives of the old 
variables with respect to the new, they involve certain combinations of 
second derivatives, but no derivatives of higher order than the second. 
Consequently the resulting invariants are not in general tensors. They 
form, nevertheless, a basis for stating an equivalence theorem. The 
sequence of invariants is composed of two invariants and their successive 
conformal derivatives, the conformal derivative being defined as the or- 
dinary derivative evaluated in a certain conformal codrdinate system. 

2. From the results of the former paper the conformal connection is 
given by 

Kh = {4} — GA) +84) - 2 enlb)/m QD 
and its law of transformation is 
Uj, = Kj us — Ki, uu + ujde + eds — gine? uit. (2.2) 
In these formulas uj = dx! /dx7, Wp = d%x'/dx dx", A (x, x) = | Ox/dx , 
¥ = log A/n, y; = dy/dx’. 

If we differentiate equations (2.2) with respect to x’, eliminate uj, by 
means of (2.2), interchange k and /, and subtract, we find as in the former 
paper equations equivalent to 


Fig: = Feo, Ui, ud up uy — oe cy + 81 Cie — Bit Ce + Bik C} + 
(5; gik — Sb gi) Vp ¥’, (2.3) 
where 


i, = Xie _ 2K 


Sa Sat + Kh, Ki, — Kf Kin, 


Cik = Vik — Vy Ve — Vp Khe, 
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and c = gc, W = ge Vik = OY;/dx*, Uj = dx'/dx’. Contrac- 
tion for 7 and / in (2.3) gives 


Fy = Fryq uf uf + (m — dew + get (nm — legjdoW, (2.4) 
where F,, = Fyqs, C = g°" Cy. Multiplication of (2.4) by g* gives 
F = F+22%n—1)c+ n(n — 1)yyV, (2.5) 
where F = pg” F,,. From (2.4) and (2.5) we find 
Vir = fir — Sq UP uk — See Veto t Khby t+ Vite (2.6) 
Fix = [Fi — 85x F/2(n — 1)]/(m —2). 
Elimination of ¥,;, from equations (2.3) by means of (2.6) leads to the 


equations, 
Chet Up = Cogs Up uh ui, (2.7) 
Chua * “el = 3 bi fit ys i fi + Sift 35 Sir fir (2.8) 
in which fi = g'? Sok Hence the first invariant of the set happens to be 
atensor. It is, except for the factor (2—m), the conformal curvature ten- 


sor of Weyl.? 
The first invariant arising from the integrability conditions of (2.6) is 


Fir = Ofjx-/ ox? — Of / dx* + Ki fon — Kh fow 


and its law of transformation is 
Sint so Sogs uf uf uj 7. Vp Chu. (2.9) 


We shall call f;,; the conformal covariant. From equations (2.9) follows: 

THEOREM 1. The conformal covariant is a tensor for spaces whose con- 
formal curvature tensor vanishes, and only for such spaces. 

If we differentiate (2.7) and (2.9) successively and eliminate each time 
the second derivatives which occur explicitly by means of (2.2) and (2.6), 
we get the desired set of conformal invariants together with their laws of 
transformation. The latter involve no derivatives of the coérdinates x 
except dx/Ox and dy/dx. 

The first invariant arising from the differentiation of (2.7) is 
Chap = OCju/Ox? + Chr Koy — Com Kfp — Cir Kby — Cieq Kfp» (2-10) 
that is, it is an invariant formed from C i in the same way as its ordinary 
covariant derivative, but with the K’s replacing the Christoffel symbols. 
This invariant is not a tensor since its law of transformation involves the 
Vi- 

Let us denote the formal expression for the covariant derivative of any 
quantity with respect to the K’s by a notation like that in (2.10). Since 
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equations (2.5) of the former paper show that (g, 2""),, = 0, we find 
from (2.8) 
Chup = Finyy + de fin/p — i Sierp + 831 8 foesy — Sin 8 Suro. 
Also by the use of the relations 
Fiujp + Fig + Finn = 0, 
Sux = Siise — firsjs 
we get 
Chup + Cire + Goan = Sefip + Si fine + 55 fia + 2 (ei forp + 


Lin far + Sirfapi) 
whence by contraction 


Chet/p = (m — 3) fiat (2.11) 


3. A system of conformal codrdinates y associated with a given co- 
ordinate system x and a point % is defined by the equations 


x! = xi, oa y' — + (Kix)o y e_ s (Kixi)o y’ y* y, (3.1) 
where 





i OK} i i j 
ja = 3 P ( oF — 2Kh, Ky + 28; fie — Sie fi), 


the symbol P denoting the sum of the terms obtained from those in the 
parentheses by cyclic permutation of the subscripts. 

We define the conformal derivative of an invariant whose transformation 
law involves only 0x/0x and dy/dx as the ordinary derivative evaluated 
in the system of coérdinates (3.1) at the origin. This derivative will have 
only 0x/dx and dy/dx in its law of transformation, that is, it is of the same 
class as the invariant which was differentiated. 

The set of conformal invariants obtained in §2 are the successive con- 
formal derivatives of Chr and fj. To prove this for a given one of the 
invariants, we differentiate the law of transformation of the invariant 
immediately preceding it in the sequence, then regard x as the coédrdinates 
of (3.1) and evaluate at the origin where the following conditions hold: 


i ee «2 
RS. fo 8 


fij = 0. 


. 


= —Ky, ¥i=0, Ki =0,— 


By a complete set of invariants we mean one in terms of which an equiv- 
alence theorem can be stated. In order to make the above set of invar- 
1 


iants complete, we add to it the conformal relative tensor*® g;;/ g” where g 
= | Li |. The statement and proof of the conformal equivalence theorem 
are omitted here because they are almost identical with those of the cor- 
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responding projective theorem. The latter is stated and proved in an 
article on ‘‘Projective Invariants of Affine Geometry of Paths,’”’ by Pro- 
fessor Veblen and the writer in volume 27 (1926) of the Annals of Mathe- 
matics. 
If we admit that such a theorem can be proved, we have the following. 
THEOREM 2. A complete set of conformal invariants for a Riemann space 
' 1 


consists of the relative tensor g;;/g", the conformal curvature tensor and the 
conformal covariant together with their successive conformal derivatives. 

If n > 3, it follows from (2.11) that the conformal covariant can be 
omitted from the list. For = 3, the conformal curvature tensor vanishes 
identically and can be omitted. 

4. By a conformally flat space is meant one which can be represented 
conformally on a euclidean space. A necessary condition for such a space 
is the existence of a preferred codérdinate system in which the conformal 
connection vanishes identically, this codérdinate system being cartesian 
for the euclidean space. The condition is also sufficient. From (2.1) we 
get on the assumption that Kj, = 0, 


m\in} = Oli} + oly} — 2 ginlin}, 


n[ jk] = gil} + gel} — gielit- 


We readily find from the last equations 


or 


Ologg; 1 O0logg 


’ 
ox* n ox* 








1 
whence g;; = g"aj;;, the a’s being constants. We therefore have g;; dx’ dx’ 
1 


= g*a,;dx' dx’, equations which show that the space can be represented 
conformally on a euclidean space. 

Since any surface is conformally flat, we assume from now on that 2 3. 

Suppose that there exists a preferred codrdinate system x such that K}, 
=0. ‘The invariants f and C vanish in such a coérdinate system and hence, 
see (2.7) and (2.9), in any coérdinate system. Conversely, if f and C 
vanish in some coérdinate system x, then equations (2.2) and (2.6) in which 
Kj, have been put equal to zero are completely integrable and define a 
preferred codrdinate system x in which K vanishes identically. 

From (2.11) it is seen that the vanishing of C for > 3 implies that of f, 
whereas for = 3, C is identically zero and the condition for flatness is the 
vanishing of f. This was first noticed by J. A. Schouten in somewhat 
different form.‘ Schouten’s conditions are 


Cirt = 0, (4.1) 


Liz: — Liz = 9, 
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where in general 
Lins — Lite = (n — 2)fia + (n — 2) 5} Chu/n. (4.2) 


My conditions are therefore equivalent to Schouten’s. It is to be noticed, 
however, that the former are expressed in terms of invariants which are 
obviously conformal in character. Schouten’s invariant (4.2) is conformal 
only when taken in conjunction with (4.1) and even then its conformal char- 
acter is not obvious. 

We now state in terms of the invariants of the present paper two theorems 
about conformally flat spaces. The first embodies a theorem discovered 
by Cotton,® the second was discovered by Schouten. 

THEOREM 3. For a Riemann space of three dimensions the conformal 
covariant 1s a tensor and its vanishing is a necessary and sufficient condition 
that the space be conformally flat. 

THEOREM 4. A necessary and sufficient condition that a Riemann space 
of more than three dimensions be conformally flat is the vanishing of the 
conformal curvature tensor. ° 

* NaTIoNaL RESEARCH FELLOW IN MATHEMATICS. 

1 These PROCEEDINGS, 11 (1925), 257-9. 

2H. Weyl, Math. Zeit., 2 (1918), 404. 

3 That this relative tensor plays a leading réle in the conformal geometry was first 
pointed out by T. Y. Thomas, these ProcEEDINGs, 11 (1925), 722-5. 

4 Math. Zeit., 11 (1921), 83; der Ricci-Kalkil, Berlin, 1924, p. 170. 

5 E. Cotton, Ann. Fac. Sci. Toulouse, 1 (1899), 412. 


GENERALIZATION OF KRONECKER’S RELATION AMONG THE 
MINORS OF A SYMMETRIC DETERMINANT 


By HELEN BARTON 
Jouns Hopkins UNIVERSITY 
Communicated March 30, 1926 
Since 1882, when Kronecker stated, without proof, that there exists a 
certain linear relation among the minors of a symmetric determinant, 
there have been many and varied proofs. It is the purpose of this paper 
to show that this relation is but a special case of a more general relation.' 


The theorem as given by Kronecker is: That among the minors of order 
m of a symmetric determinant of order 2m there exists the following linear 


relation 
| ag, | ws Dr | au | 
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where g = 1, 2,....... m h=m+1,m+2,....... 2m 
#=1,2,...m—-1l,r M=m+1,m4+2,...r—l1,m,r+1,...2m 
r=m+1,m+2,......... 2m. 


Case I. To prove this relation, let us consider the following expression, 
which involves the Generalized Kronecker Symbol.? 


Buti ..-- Bom 


m m+1 TON rong CP 928m +42 rep Im—-1 Bom—1 2g, Bom* 


For any given set of values for Bm41....8om-1; (for example, Dm, Dn+1, 

» Dm\ ar.» -Om\ 45. - -b2m—which means a complete set of b’s exclusive of 
the 2b’s which are crossed out) 8,, and Bo, must each take either of the 
values b,,4, Or bm +s. Therefore we get two expressions, namely, 


Om+r bm+1 cece bom bm +s a a b a 
oo Saree 2m i m—1 “2m—1 “bm+r om+s 
e 
bn+s 0 b 
+s "mtl ---- Yate 
~ | m mtl.... lm | oett Om—1 P2m—t Vm +s ’m+r 


which are equal in value, but opposite in sign; hence their sum equals 
zero. Since this is true for all combinations by, +1,.. .Dm\47---Om\+5---Dom 
the total expression will consist in a pairing off of terms, each pair reducing 
to zero, and thereby causing the whole expression to equal zero. 

But the expression may be written as a summation in a slightly different 
form. Consider §,, as it assumes in succession the m + 1 values m, m + 1, 
....2m; then corresponding to each value of Bm, Bm+1-....B82m Will assume 
all possible arrangements of the remaining values. But each of these ex- 
pressions is a determinant of the mth order. Therefore, if we take into 
account also the sign factor in each case, as given by the Generalized 
Kronecker Symbol, the truth of the following is evident. i 


~ Bm+1 sere ~ a, a a me 

mm+1.... 2m Bm+1~2Bm+2°°*** Bm Bom 

| ay, m+142,m+2----- am, 2m | — | ay, m 22,m+2--- m1, 2m—1 Im+1, 2m | — 
SSS AR SREP, EA Seber ken re aed eet 15 
Si | a, Le Be ee ee eee Om—1, 2m—1 Lom, a] = 0, 


i.e., | der | - pe | diz |. 


We now proceed to a slightly more general case, where instead of one 
summation symbol, occurring in the first set of values (i.e., from 1 to m) 
we have two. 














i eh! 
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Case II. We now consider the expression 


o Bm Bm+1 rie + po 


m—1l mm+1 ...: 2m 
"16m +1 [28mi2 °° °° Gm —2 Bom —2 "Bm —1 Bam—1 Bm Bom" 


This expression reduces to zero, since as before ag g,,, causes a pairing off 


of terms, so that the terms of each pair cancel one another. 

In the preceding case we arranged the expression as a summation of 
determinants, obtained by allowing £,, to assume in succession all values 
from m to 2m; so here we will group according to the values assumed by 
Bm—, and B,. For each set of values assumed by £,,_, and 6,, the remaining 
6’s will take all possible arrangements of the remaining values from (m—1, 
m....2m) and as before we shall have a summation of various deter- 
minants. But for any one set of values assumed by f,,_; and By; —Om+, 
and b,,,;, for example; there are two possible arrangements, depending 
upon whether B,,_; assumes the value D,,,, or b,4;. Since a rearrange- 
ment of values assigned to 6,,_,; and 8,, means not only a change of sign 
in the sign factor but also a change of sign of each term, it follows that the 
total result of such an interchange leaves the expression unchanged. 
Hence, there will be a complete duplication of determinants, and we may 
remove the factor 2!, provided we note that in our summation we consider 
only combinations of values for 6,,_; and 8,, and not permutations. 

In grouping these various determinants which result, we note that they 
fall into the following groups: 

1. Bm—1 and B,, assume the two values m—1 and m in the first set. 
But this is merely one combination of values, hence we represent the de- 
terminant which results by | a,, |. 

2. B»,—1 assumes values in the first set and £,, in the second and vice 
versa. In this case, suppose 8,,_; assumes the value m — 1; then let 6,, 
assume all values from m + 1 to 2m. This is nothing but case I and the 
sum of the determinants so obtained = | a, |. 

In an exactly similar way if 8,, remains fixed, with value m and we let 
Bm—1 assume all the remaining values, we again get a sum of determinants 
which equals | Op |. 

3. Both £,,_; and B,, assume values in the second set. The sign here 
is positive since there are two interchanges of values from the first to second 
set. 

Therefore, we have 


foe Bm ecee Bom 


m-1m.... | F1p4 41 22Bm4o °°: : Im —2 Bom —2 78m —1 B2m—1 78m Bom 


= [aga] — Dor lou] — Dor lanl +2 leon! = foen] — aga! — lanl +2 lapel =O 
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where g = 1,2,...... m h=m+1,m+4+2,......... 2m 
1 =1,2,....m—2,7,m j =m+1,...r—1,m—1,74+1,....2m 
¢ = 12,...m—1,7r k =m+1,...r—1,m,rt+1,....... 2m 
dike ae er ae 11, Te n =m-+1,..n—1,m—-1,14+1,..... 
to—1, Mm, fe + 1..2m 
nfs =m + i,m + 2,........ 2m. 


Hence |a,,| a pe |ap,!- 
Case III. In an exactly similar manner we see that the expression 


A= beg Bu —1 Bm cree = 


m—-2m—-l1m .... 2m 


T1641 72Bmio ++ °° Om—3 Bom —3 Bm —2 bom—2 28m —1 b2m—1 28m bom = 9- 


As Bm—2, Bm—1 and B,, assume all possible values from (m—2....2m) 
and for each set of these values, the remaining §’s assume all possible ar- 
rangements of the remaining fixed values, we get a summation of deter- 
minants. Some of these may be grouped according to case I (and there 
will be *C,, such groups); others under case II (*C2 such groups here) and 
finally a summation of determinants, in which 8,,_2, B,—-; and 6,, all three 
assume values out of the second set. 

Making use of the results obtained under case I and case II the above 
expression (A) may be written 


A= | ga | — 3|a,,| + 3 | ag, | — Do, | aus | = 0. 


Hence | apr | = x. | dy» | 


where u = 1,...... m—3, 1, 2, 3 
v =m + i1,....n—1l, m—2, 1+ 1,...%—-1, m—1, 7 + 1,... 
13-1, m, r3 + 1,...2m 
fi, fo, 3 = m+ 1, m + 2,..... 2m. 


By comparing the results of cases I, II, III, a statement for the general 
case, where there are any number of summation symbols in the first set, 
may be made. 

We note that the summations of determinants may be grouped so as to 
get + k|a,,| +>, |@yp|; the signs alternating, and k, the coefficients 
being those of the binomial expansion. Hence, we have for the case of 
four summation symbols in the first set 


|ag,| — 4 ag] + 6 |o—| — 4 )ag| + 2 low] = 0. 


Hence | aen | = a | Quy | 
where u = 1,....... 1, To, 3, 14 
> = m+ 1.05: n—1, m—3, i+ 1,....%—-1, m—2, re+1,.... 
tz—1, m—1, r3 + 1,...%4—1, m,ry + 1,....2m 
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Therefore, if we have s summations symbols in the Ist set of values, the 
general theorem is 
| ag, | am oe | ay, | 


where g = 1,2,...... m 
h=m+1,m-+ 2,....... 2m 
6 MLZ ees MS, Ty 195 005 i Wye 
v =m+i1,m+2,....n—-1,m—-s+1,n+1,....%—-1,m—-s + 


2,...-%s—1, m, 7, + 1,..2m 
ots, eee f,=m+1,m+2,......... 2m. 


The number of determinants in the summation on the right equals ™C,. 
Hence we see that Kronecker’s Theorem is a special case of the theorem 
stated here—the case where s = 1. 

1 We wish to express our indebtedness to Professor F. D. Murnaghan for suggesting 
the possibility of.extending this relation, after reading Professor E. B. Stouffer’s paper 
on this subject in these PROCEEDINGS, January, 1926. 

2 Murnaghan, F. D., Amer. Math. Monthly, 32, May, 1925 (233-241). Murnaghan, 
F. D., Bull. of Amer. Math. Soc., 31, July, 1925 (823-329). 


ON THE CLASSICAL DIRICHLET PROBLEM FOR GENERAL 
DOMAINS 


By OLIVER D. KELLOGG 
DIVISION OF MATHEMATICS, HARVARD COLLEGE 


Communicated May 17, 1926 


1. Introduction—We shall be concerned with the most general do- 
mains, or open continua, in space of three dimensions. It will be con- 
venient to denote domains, and points in them, by capital letters, and to 
denote their boundaries, and points of these boundaries, by corresponding 
small letters. We shall consider explicitly only domains, 7’, whose boun- 
daries, ¢, are bounded point-sets, a restriction which can be removed by 
inversions. 

A function F(p), defined on t, is said to be continuous if to every e>0, 
there corresponds a 6 > 0, such that |F (p) —F (q)| < ¢ for any two points 
of t whose distance, pg, < 5. By the classical Dirichlet problem, we mean 
that of determining a function, U(P), harmonic in 7, such that U(P)—> 
F(p)as P—>p. If T is not bounded, the additional restriction is placed 
on U(P) that it shall vanish at infinity like the potential of bounded 
charges. 

The problem is not always possible. Zaremba' and Lebesgue? have 
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exhibited cases in which exceptional boundary points exist, at which 
functions having otherwise the properties of solutions cannot approach 
the assigned boundary values. But in any case a function U(P), bounded 
and harmonic in 7, can be associated with given continuous boundary 
values, F(p), which will approach F(p) at every regular (not exceptional) 
boundary point.* Wiener‘ derived important properties of this function 
U(P), which we shall call the sequence solution of the Dirichlet problem. 

We shall establish the results of Wiener in a way simpler than has here- 
tofore been given, taking as point of departure a generalized function of 
Green for T, and shall obtain additional light on the sequence solution. 
We shall then give certain theorems on the capacity of point-sets, and show, 
in particular, that sets of 0 capacity are the most general at which har- 
monic functions can have removable singularities. Finally, we show that 
if the exceptional points of ¢ form a set of 0 capacity, the sequence solution 
is the only function which is bounded and harmonic in T and which ap- 
proaches F(p) at all regular boundary points.® 

2. Green’s Function.—Let [T,,] denote a nested set of normal domains 
with T as limit; i.e., a set such that the classical Dirichlet problem is 
possible for each 7;,, and any continuous boundary values, such that 7, 
is contained in all the following domains of the set, and such that each point 
of T is contained in some 7,. Then, by the reasoning of Harnack,® the 
functions of Green, G,(P,Q), for the domains, 7,, and fixed pole, Q, form 
a sequence which approaches a limit uniformly in any closed sub-region of 
T omitting Q. This limit, G(P,Q), we shall call Green’s function for T, 
leaving open the question as to whether it vanishes at all boundary points 
or not. 

From the manner of its formation, G(P,Q) is seen to have the following 
properties. G(P,Q)—1/PQ is harmonic in 7, save for a possible removable 
singularity at Q; if T is not bounded, G(P,Q) vanishes regularly at in- 
finity; G(Q,P) = G(P,Q). 

THEOREM I. G(P,Q) is independent of the particular set of nested do- 
mains, [T,,], used in its definition. For, if [T;,] denote a second set, 
G(P,Q) dominates all the functions of Green for the domains 7;,, and hence 
their limit, G’(P, Q), i.e.,G’(P,Q) S$ G(P,Q). Similarly, G(P,Q) s G’(P,Q), 
and the theorem follows. 

If T isa normal domain, G(P,Q) —> 0 as P —> }, p being any boundary 
point. The converse was first proved by Lichtenstein.* If T is not nor- 
mal, there will be boundary points at which G(P,Q) does not approach 0. 
We define as regular and exceptional boundary points, those at which 
G(P,Q) does, or does not, approach 0, respectively. 

THEorEM II. The definition of regular and exceptional points is inde- 
pendent of the position of the pole, Q. Let Q and Q’ be two points of T, 
and let R be a region consisting of a domain containing Q and Q’, together 
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with its boundary, 7, and lying in T. Let m be the positive minimum of 
G(P,Q’) onr, and M the maximum of 1/PQ on r. As G,(P,Q) s 1/PQ, 
it follows that G,(P,Q) S (M/m)G(P,Q’) on 1, and since this inequality 
holds on #,, the boundary of T,,, it holds also throughout 7,—R. Hence 
the limiting form of this inequality, G(P,Q) Ss (M/m)G(P,Q’) holds in 
T—R, and G(P,Q) must approach 0 at any boundary point at which 
G(P,Q’) does. 

The above reasoning, with the inequality given, leads at once to 

LemMA I. Let Q be confined to a closed sub-region, R, of T. Then 

(a) as P approaches a regular boundary point, G(P,Q) —> 0 uniformly 
as to Q, and 

(b) if p 1s exceptional, but such that for fixed Q’ a sequence [P,| with p 
as limit point exists for which G(P;,,Q’) —> 0, G(P;,Q) —~> 0, untformly as 
to Q,as P —> p." 

3. The Sequence Solution ——The sequence solution is derived from any 
given continuous boundary values, F(p), just as was Green’s function for 
the boundary values 1/PQ. To F(p) there correspond infinitely many 
functions, continuous in T + ¢, and coinciding with F(p) ont. Let F(P) 
be such a function. ‘Then, to each of an infinite sequence of nested normal 
domains, [T,,], there corresponds a function, u,(P), harmonic in T,, 
and coinciding with F(P) on ¢,. Wiener’s theorem is that the sequence 
[u,(P)] approaches a limit, U(P), uniformly in any closed sub-region of 
T, and that U(P) is independent of the particular extension of F(p) to the 
points of T, and of the set [T,,] employed. In calling U(P) the sequence so- 
lution of the Dirichlet problem for T and F(p), we are, of course, using 
the word solution in an extended sense, for a strict solution is not always 
possible. 

We turn now to a study of the sequence solution, starting with the case 
in which F(p) are the boundary values of a polynomial, F(P), in the car- 
tesian coérdinates of the point P of T. We also assume first that T is 
bounded. The Laplacian of F(P) is again a polynomial, 


V°F(P) = f(P). (1) 
We form the function 


, 1 
¢(P) - ae (P), bn P) =F # i I. S FQGEP,QDdVo. (2) 


In case 7}, fails to have Jordan content, the integral above is to be under- 
stood as the limit of the integral over a set of nested domains with content, 
and with 7, as limit. As f(P) is a polynominal, and as G(P,Q) S 1/PQ, 
the sequence [¢,(P)] is seen to be uniformly convergent in any closed sub- 
region of 7, and its limit, ¢(P) to be independent of the particular set 
[T,,] used. 
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THeoremM III. The function U(P) = F(P) + @(P) is bounded and 
harmonic in T, and ane the boundary values F(p) at every regular 
boundary point. 

U(P) is obviously bounded. If, in the integral for ¢,(P), G(P,Q) be 
written 1/PQ + g(P,Q), the Laplacian of ‘the first integral thus arising is 
well known to be —42f(P), while that of the second is 0. Hence V *¢,(P) 
= —f(P), and so V*¢(P) = —f(P) because the sequence [¢,4 _(P)— 
¢,(P)| is, for fixed m, uniformly convergent in T,,, and its terms are har- 
monicin7,. Thus, asV*¢(P) = —f(P),andV?F(P) = f(P),V2U(P) = 0. 

It remains to show that U(P) —> F(p), that is, that ¢(P) —>0, as 
P —> p,aregular boundary point. Since ¢(P) is independent of [7,,], we 
may assume the boundary of each 7;, to be interior to all following do- 
mains of the set. Then, for any fixed m, and any e > 0, there is a 6 > 0, 
such that when Pp < 6, G(P,Q) < «for allQ in T,. Hence, by lemma I, 
if M is the upper limit of f(P) in T, and if B is the volume of a cube con- 
taining T, 


|¢, (P) |< = ——e (3) 


for Pp < 6. On the other hand, for positive integral m, 


lenim(P) - (PSE SSS CPV SE SSS BO dVo, 


nf Sorters Tnim—Tnh 


and the last integral, extended over a region of given volume is least when 
the region is a sphere with P as center. Hence 


lenrn(P) ~ o4(P)| 3 M (E)", (4) 


where V is the inner content of T less the volume of T,. For large n, 
V, and therefore | Onim(P) — o,(P) | , is arbitrarily small, independently 
of mand P. Then, with » fixed, the right-hand member of (3) may be 
made arbitrarily small by sufficiently restricting Pp, and hence the same is 
true for ¢,4(P), independently of m. The same is, therefore, true of ¢(P), 
that is, ¢(P) —> 0 as P —> pp. 

If » is an exceptional boundary point, but one for which a sequence 
|P,], as described in lemma I (b), exists, then by the same reasoning, 
¢(P,) —> OasP,—>p. Thus the sequence solution approaches the given 
boundary values on certain sequences of points even at exceptional points, 
provided these sequences of points exist. And they do exist, save at cer- 
tain boundary points at which possible discontinuities of the solution are 
removable (p. 404). 

We have assumed, in the proof of theorem III, that T was bounded. 
If this is not the case, but ¢ is bounded, say, by a sphere with center O and 
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radius R, we proceed as follows, still assuming that F(p)-are the boundary | 
values of a polynomial, Fi(P). We form a function, H(P), as follows: 
when r = OP < R, H(P) = 1; whenr 2 2R, H(P) = 0,andwhenR <r : 
2R, H(P) is a polynomial of degree 5 in r, so chosen that H(P) and its first 
two derivatives with respect to r are continuous. The continuous ex- 

tension of F(p) to the points of T is then F(P) = F\(P)H(P). Then 
V°F(P) is bounded, and identically 0 forr 2 2R. ‘The integrals are then | 
again taken over bounded domains. 

THEOREM IV. The function U(P) of the preceding theorem is the sequence 
solution of the Dirichlet problem for T and F(p). ‘This will follow at once if 
it is clear that G,(P,Q) may replace G(P,Q) in the formula (2) without 
changing the limit, for then ¢,(P) + F(P) is exactly the sequence function 
u,(P). ‘The replacement in question, however, does not affect the limit, 
because of the uniform convergence of [G,(P,Q)] to G(P,Q) in any closed 
sub-region of T, and of the fact that G,(P,Q) $ G(P,Q) $ 1/PQ. Thusthe 
sequence solution exists, and is independent of [T,,], and approaches the 
given boundary values in the ways described above, all on the assumption 
of polynomial boundary values—a restriction which we now remove. 

Let F(p) be any continuous function, and let F(P) be a function, con- 
tinuous in T + ¢, and coinciding with F(p) on ¢. Let F’(P) bea poly- 
nomial which differs from F(P) in T + ¢ by less than e, and let u,(P) | 
and u,(P) be the sequence functions for the same set [T,] and with the 
values F(P) and F’(P) on ¢t,. Then u,(P) — ¢ S u,(P) S u,(P) + « 
in T,. As [u,(P)] approaches a limit, U'(P), uniformly in any closed 
sub-region, R, of T, for all great enough , and all P in R, U’(P)—2e S 
u,(P) Ss U'(P) + 2e. And as there is a relation of this form for every 
positive ¢, [u,(P)] approaches a limit, U(P), uniformly in R. Further- 
more, the limiting forms of the last inequalities show that the upper and 
lower limits of U(P) at each regular boundary point lie between F(p) —3e 
and F(p) + 3¢. But as these limits are independent of ¢ they coincide 
and U(p) —>F(p). The situation is similar in the case of sequence 
approach at certain exceptional points. 

Finally, the sequence solution is independent of the continuous extension : 
of F(p) to the points of T. If F(P) and K(P) are two such extensions, and 
[u,(P)] and [v,(P)] the corresponding sequences, the sequence [u,(P) — 
V,(P)] whose terms coincide with F(P)—K(P) on t, approaches 0 uni- 
formly, since F(P)—K(P) approaches 0 uniformly at the points of ¢. 
Hence the sequence solutions U(P) and V(P) corresponding to F(P) and 
K(P) coincide. 

If, in the above use of approximating polynomials, T is not bounded, we 
may take any continuous extension, F(P), of the given boundary values, 
and, using the function H(P), previously established, approximate to 
F(P)H(P) by a polynomial in the portion of T within a sphere about O 
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of radius 2R. Instead of this polynomial, we then use its product by H(P) 
for r < 2R and 0 for r 2 2R. 

4. The Capacity of Point-Sets—It is possible, as Wiener® has shown, 
to extend the notion of the capacity of a conductor to any bounded point- 
set. We shall give his definition, in slightly modified form, and give some 
properties and applications of capacity. 

Let B be any bounded point-set, and let B’ be the set consisting of B 
and its limit points. 3B’ may bound a number of domains, bit a part of 
B’, which we denote by #, will bound a domain, T, extending to infinity. 
Consider the sequence solution, v(P), of the Dirichlet problem for T and 
the boundary values 1 on ¢t. We call v(P) the conductor potential of the 
set B. The capacity of B is then the total charge of which v(P) is the po- 


: ; 1 ov 
tential, or the Gauss integral, < ps f i dS, extended over any smooth 


surface containing B in its interior, m being the normal to this surface, 
directed inward. 

THEOREM V. If Aisa part of B, the capacity of A is less than, or equal to, 
that of B. ‘The continuous extension of the boundary values may be taken 
as 1 itself, and the nested regions, [S,] and [7,,], used in constructing the 
conductor potentials w(P) and v(P) of A and B, respectively, may be so 
chosen that S, is always a part of 7,. Then, on ?#,, v,(P) = 1, while 
w,(P) S 1. Hence v(P)—w/(P) is never negative. If v(P)—w(P) van- 
ishes at any point of 7, it vanishes identically, and the capacities of A 
and B are then equal. Otherwise, let RK denote an equipotential surface 
v(P)—w(P) = C, C being small enough to insure this surface being the 


smooth boundary of a simply connected domain. On R, = [v(P) —w(P)] 


cannot be negative, since if it were, v(P)—w(P) would exceed C at some 
point outside R. Hence the Gauss integral over RK cannot be negative, 
and the capacity of A is, therefore, less than, or equal to, that of B. 

THeEorEM VI. Let T denote the infinite domain bounded by a part of 
B and its limit points, and [T,,| a nested set of normal domains with T as 
limit. Then the capacities of the boundaries, t,; of T,, converge, decreasing 
monotonely, to the capacity of B. This follows from theorem V, and the 
uniform convergence in any closed sub-region of T of the sequence [v,(P) ] 
defining the conductor potential of B, together with the uniform con- 
vergence of the sequence of normal derivatives of v,(P) to the correspond- 
ing derivatives of v(P). 

Lemma II. If Bis a bounded set of capacity c>0, the value of the conduc- 
tor potential of B at a point, P, of T, is subject to the inequalities 
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where p’ and p"’ are the points of B’ (B and its limit points) farthest from, 
and nearest to P, respectively. Because of theorem VI, and the independence 
of the sequence solution, v(P), of the nested set, [T,,], it will suffice if the 
lemma is proved on the assumption that B consists of a finite number of 
surfaces with continuous principal curvatures. Then v(P) is the potential 
of a continuous distribution of a charge, c, on B, and v(P) is decreased by 
concentrating the charge at the point of B farthest from P, and increased 
by concentrating it at the nearest point. 

Coro.tuary. If B is a bounded set of capacity 0, its conductor potential 
vanishes identically in T. For B is enclosable in the boundary, ¢,, of a 
normal domain, T,, the capacity of ¢, being arbitrarily small. Thus the 
conductor potential of B is dominated at P by an arbitrarily small number, 
that is, it is 0. 

5. Removable Singularities of Harmonic Functions ——A function which 
is bounded and harmonic in the neighborhood of a point, or of an arc of 
a smooth curve, can be so defined at such points as to be harmonic there 
also. ‘That is, such points or arcs can be the seat only of removable singu- 
larities. A complete generalization of this fact is possible by means of 
the notion of capacity. 

THEOREM VII. Let T be any domain, and let B be any part of its boun- 
dary with the following properties: (a) the set T’ = T + Bis still a domain, 
and (b) the capacity of B is 0. Then any bounded function, harmonic in 
T, may be so defined on B as to be harmonic in T’. ‘To show this, we first 
establish the following lemma. 

Lemma III. Let o be the surface of a sphere, and F(p) a function of p 
on a, bounded, and continuous except at the points of a set, s, of capacity 0. 
Then there exists a function, U(P), harmonic and bounded within o, and ap- 
proaching the boundary values F(p) at all points where this function is con- 
tinuous. We assume 0 S F(p) S$ 1. The proof can then be extended to 
the more general case by a linear transformation. We assume also that 
s is closed. Its capacity is not affected by adding its limit points. By 
theorem VI, for every positive integer n, there exists a finite number of 
smooth surfaces containing s in their interior, and of capacity less than 
1/n. Within these surfaces, a finite number of smooth curves can be drawn 
on o, bounding a domain, s,, on o, which contains s, and whose capacity 
is less than 1/n, by theorem V. Furthermore, these domains s, can be 
chosen so that each includes all the following. 

If we define F’(p) = Oons,, and F’(p) = F(p) elsewhere on a, the dis- 
continuities of F’(p) are limited to a finite number of smooth curves, so 
that Poisson’s integral defines a function u,(P), harmonic within o, and 
approaching the boundary values F’(p) save on the curves bounding s,,. 
The sequence [u,(P)] is never decreasing, and so approaches a harmonic 
limit, U(P), uniformly in any closed region within o. The bounds of 
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U(P) areclearlyOand1. Itremains to show that U(P)—>F(p) as P—>>?, 
a point of continuity of F(p). Let R be any closed region of space con- 
taining no point of s. Then, given e > O, m can be taken so large that the 
conductor potential, v,(P), of s,, is less than e in R, by theorem VI. 
But on o—Sy, Unim(P)—u,(P) = 0, for any positive integer m, and this 
difference never exceeds 1. Hence —v,(P) S uyim(P)—u,(P) S v,(P) 
on a, and hence also in its interior.1° Hence the sequence [u,(P)] con- 
verges uniformly in the closed region consisting of the points of R within 
or on og, and as the terms of the sequence have the required boundary values 
there, so does their limit. 

Returning now to the theorem, let Q be any point of B. Since T’ is a 
domain, Q is the center of some sphere, o, lying entirely in T’. Let V(P) 
be any function, bounded and harmonic in T. Then its values on o de- 
fine a function F(p) satisfying the conditions of the lemma, by theorem V. 
If U(P) be the function, harmonic in o, determined by these values F(p) 
in accordance with the lemma, V(P)—U/(P) is bounded within and on o, 
is harmonic within o except possibly at the points of B, and has the boun- 
dary values O on o except at the points of B. If the bound of the absolute 
value of V(P)—U(P) is M, this difference, and its negative, is dominated 
by the conductor potential of any set of surfaces, t,, enclosing B. Hence, 
as the capacity of B is O, V(P)—U(P) = 0 throughout the interior of o 
except on B. Hence if V(P) is defined on the portion of B within o as equal 
to U(P), it will be harmonic within «. Thus V(P) may be so defined at 
all points of B as to be harmonic in 7”. 

THEOREM VIII. Conversely, if B has the property (a) of theorem VII, 
and if any function bounded and harmonic in T can have only removable 
singularities at the points of B, then B has capacity 0. For the conductor 
potential of B is bounded and harmonic in 7, and so has only removable 
singularities. When properly defined on B, it is harmonic everywhere, 
and vanishes at infinity, and so vanishes identically. Hence the capacity 
of B is 0. 

From theorem VII, it follows that boundary sets of the type B can have 
no influence on the classical Dirichlet problem. They are called improper 
sets by Bouligand, who defines them as sets of points at which 41m G(P,Q) > 
0. It is not obvious that the concepts are equivalent, so we give a proof. 

Lemma IV. The conductor potential of any set of positive capacity has 
las upperlimit. For if its upper limit were M < 1, then, v(P) denoting this 

e conductor potential, v(P)/M would be dominated by all the terms of the 
sequence |v,(P)] whose limit defines v(P), and we should have v(P)2 
v(P)/M, a contradiction. 

Now let } be the set of points at which 4m G(P,Q) > 0. The set 7’ = 
T + 6 is an open continuum, since the points of 5 are limit points of 7, 
and since if a point of b were not an interior point of T’, it would have to 














VoL. 12, 1926 MATHEMATICS: 0. D. KELLOGG 405 


be the limit of a set of points at which 4m G(P,Q) = 0 and so could not 
itself be a point at which 4m G(P,Q) > 0. We wish to prove that b has 
capacity 0. Let A be a point of 6, and o a sphere with center at A and 
lying in 7’. Then, within and on o, G(P,Q) has a positive lower bound, 
m. Otherwise there would be in this closed region a point at which im 
G(P,Q) = 0,i., a point not in T’. Let s denote the set of points of b 
in ¢, with their limit points. Then s has capacity 0. Otherwise, the con- 
ductor potential of s would have the upper limit 1. Now on a sphere o’, 
of radius r and center Q, lying in 7, this conductor potential, v(P), has a 
maximum, M < 1, so that 3 Ga 

r\1-M 
than, or equal to 1/ron o’, and never negative. It therefore dominates the 
sequence functions G,(P,Q) in T—o’, and hence also G(P,Q). But it has 
the lower limit 0 in o, whereas that of G(P,Q) ism >0. ‘Thus the assump- 
tion that the capacity of s is positive has led to a contradiction. 

It is now possible to remove from the boundary of T a portion of }, 
contained in s but containing all the points of 5 ina sphere about A, so as still 
to leave an open continuum in place of JT. The set removed, being a part 
of s, has capacity 0. Theorem VII then shows that any function, har- 
monic and bounded in T may be so defined on this portion of b as to be 
harmonic at A. But as A was any point of b, the harmonic function may 
be defined at all points of b so as to be harmonic in JT’ = T + b. Hence 
by theorem VIII, d is of capacity 0, and so of the type of the set B. 

Conversely, if P is any point of a set of type B, im G(P,Q) > Oat P. 
For, after the definition of G(P,Q) on B has been properly made, G(P,Q) 
is harmonic at P, and so continuous and positive. 

We close with a uniqueness theorem. 

THEOREM IX. [f the exceptional boundary points of T form a set of ca- 
pacity O, there is one and only one function, U(P), bounded and harmonic 
in T, which approaches preassigned continuous boundary values, F(p), at 
every regular boundary point. One such function is the sequence solution. 
If there were two, their difference, V(P), would be bounded and harmonic 
in T, and approach O at every regular boundary point. If M is a bound 
for |V(P)|, —V(P)/M and V(P)/M are dominated by the conductor 
potential of any set including the exceptional points of ¢, and so, by a now 
familiar argument, vanish.'! 

1 Zaremba, Acta Math., 34 (1911), p. 310. 

2 Lebesgue, C. R. des Sce. de la Soc. Math. de France (1913), p. 17. 

3 Kellogg, Proc. Amer. Acad., 58, No. 14 (1923), p. 528. 

4“Certain Notions in Potential Theory,” Jour. of Math. and Phys. of the Mass, 
Inst. of Tech., 3 (1924), p. 25. 

5 Results closely parallel to a number of those to be given have recently been pub- 
lished by Bouligand in two stimulating memoirs, Mémorial des sci. math., Paris, 1926, 
and Annales de la Soc. Polonaise de Math. (1925), pp. 59-112. Of these, some had 
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been indicated to the present writer in correspondence, early in December, 1925. But 
they had been arrived at independently, in connection with his symposium on the 
Dirichlet problem at the Christmas meeting of the American Mathematical Society, 
1925, with the exception of the results on improper sets, where again, the methods 
differ. 

6 Grundlagen d. Theorie des log. Potentials, Leipzig, 1887, p. 118, ff. 

7 See Bouligand, Ann. Polonaise, p. 84. 

8 See Wiener, l. c., p. 26: P 

® See Bouligand, Ann. Polonaise, p. 104. 

1 The proof requires the fact that v,(P) approaches 1 at all points of sn. This fol- 
lows from the fact that vn(P) dominates the conductor potentials of certain spherical 
triangles inscribed sn, and the fact that spherical triangles bound normal domains. 
See Wiener, “The Dirichlet Problem,” Jour. of Math. and Phys. of the Mass. Inst. Tech., 
3 (1924), p. 104, with the addendum that parallel projections onto planes of point sets 
never have greater capacities than the sets projected. 

11 It is likely that the exceptional points of every domain form a set of capacity 0, 
in which case the above theorem would give a certain finality to the sequence solution. 
The desired result would follow from a lemma: Every closed set of positive capacity 
contains regular points. But the difficulties in the proof of this lemma have not been 
overcome. 


THE EVIDENCE FOR CHANGES IN THE RATE OF ROTATION 
OF THE EARTH AND THEIR GEOPHYSICAL CONSEQUENCES 
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This paper is a summary of one which will shortly appear in the 7rans- 
actions of the Yale Observatory and, in fact, constitutes the Introduction 
toit. The work was undertaken as a fresh attempt to see whether further 
light could be thrown on the question of the constancy of the earth’s rate 
of rotation, that is, of our measure of time. The chief evidence for changes 
of this character is derived from observations of the moon. I am not here 
mainly concerned with the secular changes due to tidal friction, but with the 
considerable fluctuations which are exhibited in the difference between the 
observed and calculated longitude of the moon when all known causes of 
variation have been eliminated. ‘The numerous investigations into ancient 
eclipses culminating in the results of Fotheringham and the work of Taylor 
and Jeffries on Tidal Friction in shallow seas are in substantial agreement 
as to the amount of the frictional effect so that it may be regarded as 
known. The tabular place of the moon due to gravitational theory is 
therefore first corrected for this effect, the fluctuations referred to being 
the differences between this corrected theory and observation. 














' 


Se a a Se ee ae ae a eh | UL 2 











VoL. 12, 1926 ASTRONOMY: E. W. BROWN 407 


Since all classes of observations of the moon in all places show these 
deviations, there is practically no doubt that they are real, that is, that 
they are not due to errors of observation, the latter being, in general, small 
compared with the magnitude of the deviations. Nevertheless, the 
differences between the results obtained from various methods of observa- 
tion and at different places, show some peculiarities which are brought out 
below for examination in order to clear the ground. 

The main issue is whether these fluctuations are due to forces, not 
hitherto recognized, acting on the moon so as to alter its motion, or whether 
they are changes in the rate of rotation of the earth. Either hypothesis 
will account for the observed data so far as the moon is concerned. Dis- 
tinction between them can only be made by observations on other bodies, 
the most important of which, on account of the material available, is the 
sun. The conclusion reached here is that the evidence favors the latter 
hypothesis. This being granted, a discussion is made of the possible 
sources of changes in the rate of rotation of the earth with such numerical 
evidence as could be obtained. 

The following brief account of the history of the discovery and of the 
attempts at explanation of these fluctuations will assist in forming an esti- 
mate of the degree of reliability of the evidence. 

Newcomb first came across the traces of the fluctuations in Part I of 
his Researches on the motion of the moon,! in which he discussed various 
observations made before the year 1750. By 1903, in a brief paper,’ he 
was able to show that they had a real existence, and six years later he gave 
a summary® of the second part of his Researches‘ which finally appeared in 
1912. He there showed that the chief fluctuation appeared to be well rep- 
resented by a harmonic term with a period of some 270 years and an am- 
plitude of about 14”. Superimposed on this were minor fluctuations of 
some three or four seconds amnlitude but irregular in period. Newcomb 
discussed the origin of these terms and, in view of the attempts to deter- 
mine the complete gravitational solution of the lunar problem, leaned 
strongly towards the hypothesis that they represented changes in the 
earth’s rate of rotation. He mentions that he tried to find confirming re- 
sults from the transits of Mercury but that he was not successful. 

The next extensive investigation was that of Cowell which was pub- 
lished in a series of papers in the Monthly Notices from 1903-5. In 
these he analyzed the Greenwich meridian observations of the moon from 
1750 to 1900. Newcomb’s results were confirmed and Cowell attempted 
an analysis of the fluctuations into two or three harmonic terms. Shortly 
after, in 1913, my theory of the gravitational problem of the moon’s 
motion being complete, and the new tables being under way, I took Co- 
well’s results and corrected them so as to make them correspond as far 
as necessary with the new theory, mainly in order to get the best possible 
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values of the principal constants. Incidentally, a comparison with New- 
comb’s results was made and a detailed determination of the minor fluctua- 
tions from both sources resulted. Within the degree of accuracy expected, 
there was good agreement in general except in the interval between 1815 
and 1845 where there were divergences which I was unable to trace down. 

In 1923, Dyson and Crommelin® took the results of Cowell and myself, 
carried them to the end of 1922, corrected them in certain particulars and 
made a full analysis of the minor fluctuations. Shortly after H. Spencer 
Jones*® published his analysis of the Cape occultations from 1880-1922, 
mainly with a view of testing the adopted constants of the theory, but also 
obtaining the run of the fluctuations during that period. A further short 
but useful series is that deduced by Russell from photographs taken at 
Harvard College Observatory during the years 1911-17.’ 

In the meantime attempts were made to see if these outstanding fluctua- 
tions could be explained as deviations of the moon’s motion. The detailed 
theory of the sun’s action on the moon and searching investigations into 
planetary and other gravitational effects made it finally clear that all such 
sources must be excluded. In 1910 I examined some hypotheses designed 
to explain them, postulating forces other than gravitational, but with 
little success. The answer to a question sometimes asked as to whether 
the Einstein effect might be responsible is clearly in the negative. One 
reason is that it is too small. A second, which applies to many similar 
hypotheses is that small disturbing forces which are functions only of the 
configuration of the system will not in general produce new periods, but 
only combinations of periods already present, and the periods which might 
be present in the fluctuations do not seem to have any such property. 

It appears therefore that we must fall back on the hypothesis that the 
fluctuations are due to changes in the earth’s rate of rotation. In 1914, 
I was able to throw some light on the question by exhibiting evidence of 
similarity between the minor fluctuations exhibited by the moon’s motion 
and the deviations of the sun and Mercury from their theoretical orbits 
during the period 1750 to 1900.8 The following year Larmor® examined 
the hypothesis that changes in the earth’s rate of rotation might be caused 
by vertical oscillations of limited areas but found that the amounts to be 
postulated were not in accordance with geophysical evidence. The same 
year, H. Glauert!® showed that a marked similarity between the fluctua- 
tions for 1865-1915 existed in the case of material derived from Greenwich 
observations of the sun, Venus and Mercury, and he also showed that their 
order of magnitude was not very different from that to be expected on the 
hypothesis of changes in the earth’s rate of rotation. In 1916, Ross"! 
showed that Mars during the nineteenth century exhibited fluctuations very 
similar to those of the sun. In 1925, Innes!* gave additional evidence from 
the transits of Mercury and the solar deviations. 
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Thus evidence has gradually accumulated which seems to justify the 
hypothesis. In view of the new material on the moon’s motion which 
has been assembled during the past few years, it seemed worthwhile to 
examine the material once again in a collected form. This first demanded 
a detailed comparison of the observations of the moon with theory, which, 
as already stated, was in large measure already available, and with this an 
intercomparison of the various series of observations to determine, if pos- 
sible, the sources of systematic differences. To test the hypothesis, the 
fluctuations should be visible in the motion of the sun, but only to one- 
thirteenth of the amount. It is necessary to combine with them the secular 
acceleration of the sun discovered by Cowell! in 1905 and later verified 
by Fotheringham" from a discussion of ancient and modern observations; 
the results were in large measure confirmed by Jeffries!®.in his calculations 
of the tidal friction in shallow seas. This demanded an examination 
of the sun’s deviations from its tabular positions and here again most of 
the material had been made. available by Newcomb and Ross. While in 
the previous attempts the effect of the minor fluctuations of the secular 
acceleration had been examined, the whole effect, namely, the total fluctua- 
tions as determined from the moon added to the sun’s secular acceleration, 
do not seem to have been applied to the observations of the sun. 

The results obtained from the present examination may be briefly sum- 
marized as follows: 

1. Observations of the sun since 1750 are consistent with the hypothesis 
that the observed fluctuations of the moon are due to variations in the 
earth’s rate of rotation, when these are combined with the known effects 
of tidal friction. In particular, the hypothesis accounts for the principal 
part of the sudden change in the deviations of the sun from its tables which 
began in 1900 and which is at present greater than 1”. The previously 
obtained correlations between these differences and those deduced from ob- 
servations of the planets are also consistent with the hypothesis but are 
partly to be explained by common systematic errors of observation. 

2. Granting the hypothesis, arguments resulting from astronomical 
data, mechanical laws and geophysical evidence point to the conclusion 
that such variations can only arise from vertical oscillations of the whole 
crust of the earth, the source of these oscillations being due to oscillatory 
changes taking place in the isostatic layer or below it. The maximum 
change of the average radius needed by the astronomical data is between 
five inches and twelve feet, according to the hypothesis made as to the 
' depth of the source. 

3. Evidence in favor of the following correlations appears to be suffi- 
ciently good for an explanation of them to be needed. 

The differences between the Greenwich meridian observations of the 
moon and the occultations as reduced by Newcomb are affected with a 
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fluctuating error which has the same principal maxima and minima as 
the minor fluctuations themselves, but the ratio of the amplitude of the 
former to that of the latter diminishes from about unity at the beginning 
of the nineteenth century to 1/8 at the present time. 

4. A fluctuation with a shorter period in these same differences appears 
to be related to the frequency of occurrence of British earthquakes. 

5. During the whole tun of the Washington observations of the sun 
(1846 to date), the differences between them and the tabular values follow 
the same major fluctuations as the differences between the Greenwich 
observations of the sun and the tables. The change appears to have a 
period of about forty years. 

6. The large deviations of the Greenwich observations of the sun 
from the tabular values between 1805 and 1825, taken in conjunction with 
the similar large deviations of the Greenwich meridian observations of the 
moon from the occultations, point to the conclusion that the meridian ob- 
servations are defective in this interval and that it is unsafe to make use 
of them with full weight for the discussion of tabular differences, unless the 
source of the apparent errors can be found and the errors, if existent, 
corrected. The indications are that for discussion of variations of the 
earth’s rate of rotation, the occultations rather than the meridian observa- 
tions of the moon should be used in this period. It follows that harmonic 
analyses based on the meridian observations before 1850 are of doubtful 
reality and that the minor fluctuations are much more irregular in character 
than had previously been suspected. 

7. It is doubtful whether attempts at harmonic analysis of the fluctua- 

tions will give any clues as to their origin. 
8. A somewhat doubtful correlation is that between the Greenwich 
(and Washington) tabular differences for the sun and the differences be- 
tween the Greenwich meridian observations of the moon and the occulta- 
tions since 1870. Another, which may be accidental, between the same 
series appears before 1840 but it requires a lag of ten years in the former in 
order to make it correspond with the latter. 

In the discussion on the forces which can change the earth’s angular 
velocity, where I have attempted to prove that the only hypothesis which 
can account for the facts is an oscillatory change in the earth’s mean radius. 
two calculations have been made. In one the earth is supposed to expand 
and contract throughout its whole mass uniformly. In the other, all the 
expansion and contraction is supposed to be due to changes of pressure in a 
layer 50 miles deep, changing the height of the crust above it. These are 
set down merely as outside limits which are in accordance with the knowl- 
edge we have at present concerning the earth. The limits can possibly be 
narrowed, but I am not concerned here with attempts to do so nor have I 
tried to make any suggestion as to the physical or chemical mechanism by 
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which such expansions or contractions can be produced. As I view the 
problem, the argument depends on a chain of reasoning which appears to 
require the existence of some such force well below the crust of the earth. 
The strength of the chain depends on the strength of its weakest link. 
The chief steps are as follows: 

1. The astronomical data. There can be little doubt that the principal 
part of the apparent variations in the moon’s motion are real and not due 
to errors of observation. 

2. The evidence that the observed variations are variations in the rate . 
of rotation of the earth and not in the rates of revolution of the moon, sun 
and planets about the earth. The greater part of the paper is devoted to 
setting forth the evidence for this, involving attempts to distinguish be- 
tween systematic errors of observation and real differences between gravi- 
tational theory and observation in the motions of these bodies, especially 
in that of the sun. 

3. The arguments by which sources for these variations are traced 
to forces well below the outer crust of the earth. The evidence is largely 
negative, that is, it consists of attempts to show that, unless matter acts 
on matter in some way at present quite unknown, the forces cannot proceed 
from any other source. 

4. Assuming the last result as provisionally established, the collateral 
evidence which can be brought forward to support it in the form of correla- 
tions between geophysical and the astronomical data, and simplified ex- 
planations of geodynamical questions. 

An oscillation with a period between 200 and 300 years in the frequency 
of Chinese earthquakes brought forward by H. H. Turner!® in 1920, and 
compared by him with the principal fluctuation in the moon’s motion 
may be significant in this connection. 

As the main argument must ultimately depend on the numerical data, 
the latter have been set out in detail as far as possible with full references 
to their sources, and with an exact account of what changes or corrections 
have been made. The treatment of these data has necessarily been such 
as to bring out as prominently as possible the facts needed for the argu- 
ment. That it also brings out other results for which explanation is not 
forthcoming is inevitable. The question as to whether the latter are rel- 
evant to the main issue is discussed in some detail, _ 

While the discussion of the effect of changes in the external radius of the 
earth would lead to the expectation of correlation between geophysical 
and astronomical data, the mechanism which is responsible for a correla- 
tion between the errors (if existent) of the Greenwich observations and the 
frequency of certain seismic data is not apparent. Owing to the lack of 
uniformity of the earth’s crust at least down to the isostatic layer of com- 
pensation, we may expect differential vertical motions of the crust super- 
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imposed on the main vertical oscillation common to the whole earth. 
The magnitude of that oscillation exhibited by the moon which is corre- 
lated below with the frequency of occurrence of British earthquakes seems 
to be rather greater than can be attributed to changes of the vertical at 
Greenwich. Such an explanation would seem to demand corresponding 
changes in the sun’s motion, and there is not much, if any, evidence of them. 
If, however, as seems likely, the effect on the latter is only '/;; of that on 
the former their absence would be explained. The well-established corre- 
lation of the differences between the Greenwich meridian observations 
and the occultations, with the differences between the observations and 
theory would seem to result from the same cause, but it is difficult to sug- 
gest any reasonable hypothesis to account for the manner in which it is 
brought about. 

Attempts were made to connect these differences with those obtained by 
Sampson!’ in the time determinations at different observatories. His 
proof that, at least at Edinburgh,'® they were connected with level errors 
seemed to hold out some promise and I made attempts at correlation be- 
tween his results and those simultaneously furnished by observations of 
the moon and sun, but obtained nothing of sufficient significance to re- 
produce here. 

A theory by Prof. Joly, reported in the Observatory for 1926, February, 
also postulates vertical oscillations of the whole crust.'? The basis appears 
to be its value in explaining the phenomena of dynamic geology. He 
believes that the cause is to be found in the thermal effects of radium act- 
ing in a substratum of basalt at a depth of some 30 kilometers and that it 
is measured in time by geological periods. This basalt was alternately 
liquid and solid; at the present time it is solid. In this paper I deduce 
from astronomical data, secured during the past two centuries, osciJ!ziions 
measured by years and by one or two centuries, and place the ovigin at 
least at a depth of 80 or 90 kilometers, but make no hypothesis as to their 
cause. The nature of the oscillation which he deduces is therefore quite 
different from that of this paper. 

1 Wash. Obs. App., 2, 1875. 

2 M. N., vol. 63, p. 316. 

3M. N., vol. 69, p. 164. 

4 Amer. Eph. Papers, vol. 9, pt. 1. 

5 M. N., vol. 83, p. 359, and App. G. of the Greenwich Observations for 1920. 

6 M. N., vol. 85, p. 21. Cape Obs. Annals, vol. 8, pt. 8, 1925. 

7A. J., vol. 36, p. 153. 

8 Brit. Assn. Report, 1914, p. 311. 

9 M. N., vol. 75, p. 211. 

0 M. N., vol. 75, p. 489. 

11 4. N., vol. 225, p. 109. 

122.4. J., vol. 29, p. 152. 

13 M. N., vol. 66, p. 3. 














VoL. 12, 1926 PHYSIOLOGY: C. LADD-FRANKLIN 413 


4M. N., vol. 81, p. 104. 

1 Trans. R. S., vol. 221A, p. 239. 

16 M. N., vol. 80, p. 807. 

WM. N., vol. 81, p. 89. 

18 M. N., vol. 85, p. 560. 

1” A summary of his book, The Surface of the Earth. 


THE REDDISH BLUE ARCS AND THE REDDISH BLUE GLOW 
OF THE RETINA: SEEING YOUR OWN NERVE CURRENTS 
THROUGH BIOLUMINESCENCE 


By CHRISTINE LADD-FRANKLIN 
CoLuMBIA UNIVERSITY 


Read before the Academy April 27, 1926 


A simple band of bright red light thrown upon a screen in a dark room 
gives rise to a very curious phenomenon—discovered in the first instance 
by Purkinje.! What you see on the screen is not only the red band but 
also, projecting out from it on both sides, big slightly reddish blue arcs. 
They are not of the color of the visual purple, which is a slightly bluish red. 
If one considers the shape and the angular size of these arcs, it is perfectly 
plain that what one sees, as an entoptic phenomenon, is certain fibres of the 
optic nerve which lie on the surface of the retina and which proceed to 
their exit point, the papilla. But why are they visible? The explanation 
that has been given of this phenomenon hitherto, by Gertz, by Troland 
and by others, is that the nerve current by which one sees the red band 
gives rise to a secondary induced current in adjoining nerve fibres. 
Such a current as this, however, would not be provided with the right 
“place coefficients’’—it would not enable one to see the stimulated fibres 
in the place where they lie. Let me put this hypothetical question: 
Suppose that one could lay open an optic nerve fibre all the way from the 
retina to the occipital lobe of the cortex, and suppose that one were to 
pinch it in half a dozen different places. The pinching of this nerve fibre 
would give one a light sensation (either chromatic or achromatic), but 
where would that sensation be seen to be ?>—not, certainly, in the half dozen 
different places in which the nerve fibre is pinched. The sensations would 
one and all seem to be in that direction in the external world which is in 
one-to-one correspondence with the rod or cone in the retina from which 
the fibre which is pinched hascome. We have the sanie situation when the 
man who has had his leg amputated at the knee has a nerve fibre pinched— 
he feels a sensation of tickling (or something), but he feels it not in his knee 
but in the ball of his foot.? It is therefore plain that in the case of the 
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Blue Arcs of the Retina no possible nerve current evoked in adjoining fibres 
by the fibre which is conveying the red-light exertation would enable you 
to ‘‘see’’ those fibres, the actual place-coefficients attached to a sensation 
produced in this way would be entirely wrong. 

There is another insuperable difficulty in resorting to any sort of a nerve- 
fibre stimulation as the cause of the reddish blue arcs. An after image of 
these blue arcs can be obtained; one gets a sensation which meets the re- 
quirements of the after image both in achromatic intensity and in chroma— 
it is of a slightly greenish yellow. But an after image does not occur after 
a stimulation of the visual mechanism by an electric current. This one 
would have been inclined instinctively to take as self-evident; but it has 
now been proved beyond question by the work of Lazareff, in which he 
finally demonstrates that nerve action is not (and cannot be) subject to 
fatigue.® 

That nerve fibre when stimulated gives off some sort of an emanation 
which enables it to take its own photograph has now been proved by No- 
dan;‘ and it is perfectly easy to suppose that this ‘“emanation’’ may be 
of the nature of visible light or at least of a nature to produce visible light 
by means of fluorescence, which is known to occur in the retina. 

Various details of this phenomenon—the blue glow, which has been en- 
tirely overlooked by most observers, the absence of the effect in the case 
of individuals who have myelinated fibres in this region of the retina, and 
other considerations—make it very certain that its real cause is a bio- 
luminescence on the part of the stimulated nerve fibres. It would follow 
that all nerves when stimulated shine by their own light—a light which is 
invisible, of course, when the nerves are not non-myelinated. One has not 
discovered this before, because, although one devotes much time to study- 
ing the stimulated nerve fibre, one does not do it in an absolutely dark 
room. 

1 This rare work of Purkinje—so rare that Gertz reproduces what Purkinje says 
about this phenomenon, because so few of his readers will be able to see his book—has 
now been reproduced in Czecho-Slovakia. Purkinje, Johann Evangelista, Opera Umnia. 
Praze: C. Calve, 1919. 

2 To express this situation the psychologist has need of a new term. It is absurd to 
apply the phrase “‘local sign’”’ to the case of the knee and the foot. ‘Sign’ is a wrong 
term anyway, for there is every reason to think that we have here to do with an actual 
element of consciousness, and not at all with a sign simply. I propose ‘‘place- 
coefficients” as the proper term for this form of attachment to the pure sensation. 

3 Lazareff, Comptes rendus de l’ Academie des Sciences, 1925. 

Katz, David: “Die Illusionen der Amputirten,” Betheft zur Zts. f. angewandten 
Psychologie, Leipzig, Barth, 1921. 

4 Nodan, Alfred: Jbid., 1923-1926. 
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THE METABOLIC EFFECT OF ENEMATA OF ALCOHOL, 
DEXTROSE AND LEVULOSEIN HUMANS 


By THORNE M. CARPENTER 


NUTRITION LABORATORY OF THE CARNEGIE INSTITUTION OF WASHINGTON, BosTON, 
Mass. - 


Communicated April 28, 1926 


The usual path of ingestion of food and drugs by man is by way of the 
mouth. Under some conditions, such as congenital or accidental ob- 
struction of the alimentary canal, unconsciousness, anesthesia or un- 
willingness of the individual to coéperate (insanity), it is necessary to re- 
sort to other paths of ingestion. Those commonly used are the intra- 
venous, subcutaneous, intraperitoneal and rectal injections of the substance 
needed or under observation. Any study of the metabolic effect of one of 
these methods of ingestion has therefore a practical significance, as it re- 
sults in an increased knowledge regarding the effectiveness of the method. 
Such an investigation has another interest. It is conceivable that when a 
substance is introduced into the circulation directly or ingested at a differ- 
ent point in the alimentary canal (e.g., the rectum), its path of distribution 
may be different than when the ingestion is by mouth, and the character 
and speed of metabolism may also be altered. 

This communication has to deal with the effects of enemata of ethyl 
alcohol, dextrose and levulose upon the metabolism of humans observed 
in an investigation which began as a part of the program! of study of the 
Nutrition Laboratory on the physiological and psychological effects of 
ethyl alcohol. Other experiments were added in order to include for com- 
parative purposes the effect of rectal injection of solutions of dextrose and 
levulose. 

The subjects were four medical students presumably in good health. 
Problems studied were how completely and rapidly the materials were 
absorbed, what effect, if any, the absorption of these materials had upon 
the volume and some constituents of the urine, whether the pulse-rate was 
changed and to what extent the material absorbed was metabolized as 
indicated by a study of the respiratory exchange... The last was the chief 
problem of the investigation. The interpretation of the results depended 
in the main upon the course of the respiratory quotient. 

The respiratory quotient is the relation by volume of the carbon dioxide 
eliminated to the oxygen absorbed when a substance is oxidized. For the 
utilization of fat the value is 0.713, for carbohydrate 1.00 and for the 
portion of protein oxidized in the body it is 0.81. The respiratory quotient 
of alcohol is lower than any of these, namely, 0.667, so that if alcohol is 
metabolized in the body in such a manner that no.other excretory products 
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are formed than carbon dioxide and water, the result will be a lower re- 
spiratory quotient than that existing before alcohol was ingested. That 
this occurs when alcohol is given by mouth has been shown by a previous 
investigation in this Laboratory.* Oxidation of dextrose or levulose, 
on the other hand, gives a respiratory quotient of 1.00 and ingestion by 
mouth results in a rise in the respiratory quotient,’ the greater and quicker 
change coming after the taking of levulose. 

The absorption of the substances used in the enemata of the present 
study was determined indirectly; that is to say, after a solution had 
been retained a given time in the rectum and large intestine, a cleansing 
enema was given to remove the unabsorbed material, and to be sure that 
this removal was complete, two and sometimes three such enemata were 
used. The amounts of the unabsorbed substance thus secured were then 
determined by suitable chemical methods. 

Urine was collected in the period before an experiment (preliminary) 
and for the time just prior to and following the rectal injection (experi- 
mental). The volume, total nitrogen and chlorides were determined for 
these urines. When alcohol was injected, the amount eliminated by the 
kidneys was also estimated. Ina number of special experiments, the urine 
was collected in short periods. Determination of the alcohol content of 
these urines was expected to throw light on the metabolism of alcohol. 
For comparison purposes there were a few experiments in which similar 
amounts and concentrations of alcohol were taken orally and the alcohol 
determined likewise in the urine collected in short periods. The pulse- 
rate was taken periodically in all experiments in which the respiratory 
exchange was measured. The latter was determined by means of the 
gasometer method, in which the expired air is separated from the inspired 
air by valves and the expired air is collected in a gasometer and analyzed. 
In several series of experiments a chamber‘ respiration apparatus was 
used. A 0.6 per cent solution of sodium chloride was used alone in a num- 
ber of experiments as a control solution to determine the effect of rectal 
injection of a liquid upon the several physiological and chemical processes. 

Ethyl alcohol is very rapidly and nearly completely absorbed when given 
inenemata. In 31 experiments in which quantities of alcohol varying from 
11.0 grams in 220 cc. of 5 per cent solution to 60.8 grams in 810 cc. 
of 7.5 solution were injected, there were only two cases in which the ab- 
sorption was not better than 98 per cent, and in these two instances the 
duration of the periods during which the solution was retained was 1 hour 
and 4 minutes, and 2 hours and 3 minutes. Alcohol was almost invariably 
found in the urines collected. ‘The maximum amount eliminated in this 
manner was only 1.9 per cent of the alcohol given. This occurred in one of 
the experiments in which urines were collected in short periods and as of- 
ten as possible, after the giving of 37.5 grams of alcohol in a 5 per cent solu- 
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tion, the collection lasting about 7 hours. The quantity of urine eliminated 
after alcohol was given was nearly always increased over that observed in 
the preliminary period before the experiment. In spite of this, the ni- 
trogen and chloride eliminated were decreased after alcohol, the nitrogen 
slightly and the chloride to a more marked degree, indicating that the 
effect of the alcohol was to promote the elimination of water without re- 
thoving other constituents of the tissues at the same time. 

The most striking effect of the rectal injection of alcohol was the marked 
rise in pulse-rate. Previous researches of Dodge and Benedict,’ Higgins? 
and Miles* have shown that the effect of alcohol is either to prevent the 
natural fall in pulse-rate, which would take place in the course of an ex- 
periment if no alcohol were given, or to produce a slight increase in the 
pulse-rate. In other words, there was usually a relative acceleration. 
This amounted to only a few beats at the most. In the experiments. with 
rectal injection, however, the increases in pulse-rate were of much greater 
magnitude. The most marked effect observed was in a series with a sub- 
ject who slept in a respiration chamber for the entire night after the rectal 
injection of 37.5 grams of alcohol in 500 cc. of water. In the five ex- 
periments of the series, the pulse-rate rose in the course of 6'/. hours from 
an average level of 62 to an average of 73, and in an individual experiment 
it rose from 66 to 87 in 61/2. hours. Rectal injection of alcohol has a more 
specific influence upon the pulse-rate than the giving of like quantities by 
mouth. 

The respiratory exchange indicated very clearly a utilization of alcohol. 
The respiratory quotient began to lower within two hours after the rectal 
injection of alcohol and the greater the quantity of alcohol given, the more 
marked and persistent was this effect on the quotient. In the experiments 
by Higgins? the decrease in the respiratory quotient commenced much 
sooner. Also in experiments in this research when alcohol was given by 
mouth, the respiratory quotient was lowered within one-half hour. The 
studies of the urine, in which the alcohol elimination was determined by 
periods, indicated that the peak of concentration occurred at about the 
same time whether ingestion was by mouth or by rectum, that is, within 
the first hour and one-half, and that the alcohol disappeared from the 
urine in about the same length of time for equal quantities injected, i.e., 
25 grams in about 5 hours, an average disappearance by various ways of 
5 grams per hour. 

Calculated from the respiratory exchange and by means of Rosamann’s 
formula,’ the maximum amount theoretically metabolized after rectal in- 
jection of alcohol was 4.4 grams per hour when 37.5 grams were given in 
500 cc. of water. There is an unexplainable inconsistency in the re- 
sults of the respiratory exchange measurements and the alcohol in urine 
‘ determinations in that the respiratory quotient in these experiments re- 
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mained much depressed at a time after the giving of alcohol when urine 
analysis would indicate a tendency for the alcohol to disappear. 

The solutions of dextrose and levulose in 5, 6 and 10 per cent concen- 
trations were not so rapidly or so completely absorbed as similar solutions 
of alcohol which were retained for the same length of time. The maximum 
absorption of dextrose was 90 per cent in the case when 30 grams were 
given and the material injected was retained about 6'/, hours. The maxi- 
mum amount absorbed per hour was 7.8 grams in a period of 2'/, hours. 

Levulose was, in general, absorbed slightly more rapidly than dextrose. 
The apparent absorptions are higher on the average than those for dex- 
trose, several of them being over 90 per cent with a 5 per cent solution. 
When 50 grams in solution were retained 4'/. hours, an average of 10.4 
grams per hour were absorbed. Indications from the results are that the 
rate of absorption of dextrose and levulose is greatest at first, and prolong- 
ing the time of retention of the material injected does not increase pro- 
portionately the amount of absorption. 

In the majority of experiments with dextrose the respiratory quotients 
were definitely raised, thus indicating an increased utilization of carbo- 
hydrate in the metabolism. The maximum increase in carbohydrate 
utilized calculated from the changes in the respiratory quotient was 4 grams 
per hour. If all of the carbohydrate utilized after the injection of dex- 
trose was supplied by the dextrose absorbed, the absorption recorded would 
hardly suffice for more than two or three hours. If one considered only the 
increase in utilization of carbohydrate, the amount absorbed would be 
sufficient for a much longer time. 

On the other hand, the rectal injection of levulose does not, in the ma- 
jority of cases, produce any marked rise in the respiratory quotient. This 
is in direct contrast to its effect when introduced by mouth. When it is 
introduced by mouth, it has a greater effect and produces a more prompt 
increase in the respiratory quotient than does an equal quantity of dex- 
trose. Levulose, however, of all of the substances given by rectum, caused 
the greatest decrease in nitrogen elimination in the urine, indicating a 
sparing action upon the protein metabolism. This was not the case with 
dextrose, which had an effect upon the nitrogen in the urine about equal 
to that due to alcohol. Dextrose produced a greater increase in the 
volume of urine than did levulose while, on the other hand, dextrose had less 
effect upon the elimination of sodium chloride. In other words, with 
dextrose there were indications that the liquid injected and sodium chloride 
were eliminated while with levulose, the liquid injected and sodium chloride 
were retained. 

The results with alcohol, as well as with dextrose and levulose, indicate 
that the metabolism of materials introduced by rectum is different than 
that produced when they are given by mouth. The usual explanations - 
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of the different effect of rectal injection upon metabolism are that the 
material is absorbed in such a manner that it goes directly into the systemic 
circulation and avoids the portal circulation and that differences observed 
in the results are due to the slowness of absorption. ‘The first explanation 
seems untenable because all of the indications both from this work and from 
other researches are that the material passes high enough into the large 
intestine to be absorbed into the branches of the portal circulation and that 
consequently it must pass through the liver by means of the portal circula- 
tion. Experiments by others with small quantities of sugar given by 
mouth indicate that it was not the small quantity absorbed which produced 
the effects in question so that the slowness of absorption can hardly ex- 
plain the results. From the results in this investigation it is believed that 
the material introduced by rectum passes into the portal circulation and 
through. the liver and then throughout the body. As a hypothesis it is 
suggested that the liver, as well as the rest of the alimentary tract, under 
the conditions of the experiments is not in an active state so that the 
substance is metabolized throughout the body in a manner similar to the 
utilization of material in muscular work. Dextrose and alcohol are both 
substances which can be utilized in the body when given in enemata. 

A review of the earlier literature, the details of the experiments, a dis- 
cussion of the practical applications and theoretical calculations, together 
with a suggested hypothesis, appear in Publication No. 369 of the Car- 
negie Institution of Washington. 


1 Benedict and Dodge: Tentative Plan for a Proposed Investigation into the Physio- 
logical Action of Ethyl Alcohol in Man. Proposed Correlative Study of the Psychological 
Effects of Alcohol on Man. Privately printed, Boston, Massachusetts, 1913. See, also, 
Dodge and Benedict, Carnegie Inst. Wash. Pub., No. 232 (1915), pp. 266-275. 
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6 Miles, Carnegie Inst. Wash. Pub., No. 333, 1924, p. 225. 
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ERRATA 
These PROCEEDINGS, pp. 335 ff., May, 1926. 
3n 


3n 
In equation (1), >, should read >. , 
j=l j=4 
3n—b 3n—b 
in equation (2), 7 should read 2 , 
j=1 j=4 


in equation (3), »© should read vf), 


and on page 337 “k the gas constant per gram” should read ‘“‘k the gas constant per 
molecule.” —E. O. SALANT 

















